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I. INTRODUCTION 



1. Purpose . Increasing prices and decreasing supplies of 
fuel oil has created a need for highly efficient ship pro- 
pulsion systems capable of utilizing fuels other than oil. 

The purpose of this thesis is to perform a preliminary design 
of a combined cycle propulsion plant equipped with a multi- 
fuels combustor capable of burning coal or oil. 



2. System Description . The system investigated is a closed 
regenerative Brayton Cycle combined with a bottoming Rankine 
cycle. Working fluids for the cycles are air and water 
respectively. The heater for the Brayton cycle is a pressur- 
ized fluidized bed (PFB). The heater cycle is essentially an 
open regenerative 3rayton cycle. This thesis was conducted 
considering the heater only in its coal burning mode of opera- 
tion. Work done by the Brayton and Rankine cycles will drive 
a common shaft through a double reduction type gear. An enthalpy - 
enthropy diagram is shown on Figure 1 . 



3* Scope . 

a. Objectives. The objectives of this thesis are to: 
(1) Perform a thermodynamic analysis of the system 
for various closed Brayton cycle pressure ratios utilizing 
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Figure 1 . 







7 



realistic values for component efficiencies and losses. 

(2) Estimate major component volumes and thus obtain 
a total volume for set of parameters listed in (1) above. 

(3) Select a set of parameters for design based 
upon the above results. 

(4) Discuss the technical feasibility and implications 
of the design as a means for naval ship propulsion. 

b. Design Criteria . 

(1) The system parameters selected for design were 
those which gave the lowest total volume within an acceptable 
range of thermal efficiency. 

(2) An acceptable range of thermal efficiency is 
defined as t 1% of the maximum efficiency. 

c . Limitations to Scope . 

(1) No component or system testing was carried out. 

(2) The system was designed with the aid of a 
computer utilizing performance parameters of currently avail- 
able equipment. 

(3) The concept of pressurized fluidized bed (PF3) 
combustion has not been proven in a shipboard environment. 

The thesis was conducted based on the assumption that a PFB 
system would be put to sea in the near future. 
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II. THERMODYNAMIC ANALYSIS 



1 . The Brayton Cycle 

a. Working cycle . This is a closed Brayton cycle 
utilizing air as the working fluid. Air was chosen since 
this is proposed for surface ship application and would thus 
eliminate the need for carrying a reserve supply of gas if 
something like Helium or CO 2 were used. 

Advantages of a closed cycle over an open cycle are 
the followings 

(1) Compressor suction pressure is no longer limited 
to ambient. By increasing the working pressure of the cycle 
a significant reduction in engine size can be achieved for a 
specified output. 

(2) Higher working pressures also increases the heat 
transfer capacity of the working fluid thus resulting in 
smaller, more effective heat exchangers. 

( 3) -he use of a high pressure fluid gives the added 
advantage of regulating the cycle pressure level by varying 
the inventory of gas in the loop. If this is done in a manner 
which kept the maximum temperature and speed constant, the 
turbomachinery vector diagrams would remain essentially un- 
changed. The result is little or no change in compressor or 
turbine efficiencies with varying load, thus maintaining high 
part load efficiency. 

A Temperature - entropy (T-s) diagram of the working 
cycle is shown in Figure 2. 
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Figure 2. 



point Is 
2s 

3: 
4s 
5s 
6 $ 
7 : 



compressor suction, cooler discharge 
compressor discharge, regen inlet (cold) 
regen outlet (cold), heater inlet 
heater outlet, turbine inlet 
turbine outlet, regen inlet (hot) 
regen outlet (hot), WHB inlet 
'.‘JHB outlet, cooler inlet 



QgS heat input from heater cycle 

t 

Q^s heat released to Rankine cycle 
Q c $ heat released to sea water cooler 
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The analysis was conducted by assuming values for 
T^, T^, compressor polytropic efficiency), (P) (turbine 

polytropic efficiency), (regenerator effectiveness), and 
(total pressure drop). 

The specific heat was assumed to vary with temperature 
according to the following relation: 

Cp(T ) = .2475 - ( 3.759x10“-’) T + (5 . 106x1 0 - 8 )T 2 -( 1 .231xlO -11 )T 3 

for T in °R 

Cp in BTU/lbm-°R 

The following are defined : 

compressor pressure ratio ( Pg/P^ ) 
turbine pressure ratio (P^/P^) 

% total pressure drop around closed Brayton loop 

power (net) from closed Brayton cycle 
mass flow of working fluid in Erayton cycle 

temp in °R at point i 

enthalpy at point i 

gas constant 
pressure abs . at point i 

specific volume 
regenerator effectiveness 

The basic procedure was to determine Tg, T^, T^, T^, 

Q^ n , and Wg for various pressure ratios ( TT* ) . 



7T„ 

v- 

TTt 
wi ' r 



k 



R 

P. 

l 

V 

E„ 
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Determination of (ref. Dixon): 

Small stage polytropic efficiency is 

, p) = ihi z zh. 

c (p; dh C P At 

since for an isentropic process = 0=dh^ s ~ vdp 

Substituting V = RT/p gives 

i )c(f 3 ) = £L 1 dp 

Cp * p * dT 



thus, 



-71 



cl' = 



T, 



7 " 



o 

/«■ 



A_ ( 

(p; y 



<^p 



in 






pc (p) l p. 



By substituting the equation for C = C (T) one gets 

P P 

.2475 ln/T 2 \ 3.759xlO" 5 (T 2 -T 1 ) + ( 2 .553xl0 -8 ) (T £ 2 -T 

\*lj 

(4.103x10“ 12 ) (T 2 3 - T^) 



= R 



Mr) 



7- /o 7T C 



T 2 can then be determined iteratively 



Determination of T> : 



In a manner similar to the compressor above, 



,T C 



■^M T = R /> 



% 
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After integrating, this becomes 



.2475 ln | T 5 j- 



(3.759x10' 5 )(t 5 -t 4 )+(2.553x10‘ 8 )(t 5 2 -t 4 2 )- 
(4.103x10~ 12 ) (T^ 3 - T 4 3 ) 



* R %(<’)(" K) 




can then be determined iteratively. 



Determination of T ~ and Tg: 

and can be found from an assumed regenerator 
effectiveness (E R ) using the following relations 
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b. Heater Cycle . This is an open Brayton cycle also using 
combustion air as the working fluid. The heater itself is 
a pressurized fluidized bed utilizing coal as fuel. A T-s 
diagram of the heater cycle is shown in Figure 3* 

Major components of the cycle consist of a heater, re- 
cuperator and turbine driven compressor. The turbine is 
powered by exhaust gases leaving the recuperator. 

The analysis was conducted by assuming values of T^» 
^c(P)» y.f(P), and 0 (equivalence ratio). Both the 

exhaust products and combustion air were assumed to be perfect 
gases with specific heats varying in accordance with the re- 
lation given in l.(a). 



TT-th 

< 7 % 

m 



The following are defined ; 

7T C h - heater compressor pressure ratio P 2h/ P lh ) 

heater turbine pressure ratio ^ P ?k //P 6h^ 

% total pressure drop around heater loop 
mass flow of air in heater cycle 

mass flow of fuel in heater cycle 

stoichiometric fuel - air ratio of fuel used 

Lower Heating Value of fuel used 
absolute temp at point i^ 

enthalpy at point i^ 

efficiency of the heater cycle 
recuperator effectiveness 



( f/a ) s 

LHV 

T ih 



h ih 

H 



V 



J REC 
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T 



ambient .. 




S 



Figure 3 ' 



point 


Ihs 

2h: 




3h: 

4h: 




5h: 

6h: 




7h: 


V 


heat 


« 

Q. : 
m 


heat 



compressor inlet 

compressor outlet, recup inlet (cold) 
recup outlet (cold), heater inlet 
max air temp in heater 
heater outlet, recup inlet (hot) 
recup outlet (hot), turbine inlet 
turbine exhaust 



heat released to working cycle 



LHV 
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Assume the heater cycle to be a control volume (C.V.) 



as shown in Figure 4. 



m f LHV 




I 

I 

I 



C.V. 



m 



a 



h 



Ih 







V 

• • 

(m a + m f ) 



h 



7'n 



Figure 4. 

Qg leaving the C.V. must be equal to Q 3 entering the 
working cycle thus, 



= ^3 ^4 " + ^f^ h 4h “ h 5h^ 



Performing a heat balance on the C.V. gives us 



+ (m a + m f ) hr, h = (LHV) + m„h 



a lh. 



Therefore, 






<*3 Q& 



Q:h fa-? CW 



-[(i +Ws 0)^ - U /[(%) s <t> LHV] 
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From a power balance between the turbine and 



compressor 



^c = W T ^ no ne_t work ) 



m 



pc 



(T 



h2 



- T hl> ■ (“a 



* V C pT (T 6h - T 7h> 



but 

Lh =T k (iTc)^^ 

T‘'" = hk fa*d W '° TT^ * (l I ■ (?%, 7 tt c k ) 

thus , 

c* - ! ] = [ 1 1 (%)s <Pj G, t 7 , * - ,] 

= [I + faA n C t , T« { [ i - (f J n ] c * ^ [(n^ * - 1 



or, 



G ^ * [ I + fads n C r , T % = 



[< 










Qv-yt- 
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therefore, 



( ' * p] Cfi T 1L (mk) ** » 

For minimum volume heat exchangers we want the largest 
f AV/ \ 

( /P/fix compatible with the above eqn. 



differentiating, 






Setting / A ? ) 

\ ?/Th 

"^ch 



0 and solving forTTch gives 




[i +{ r A)t $ ] 9r ft fa % - 

2/v T\ ( /{V - Cfn. fi) 
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By substitution 



A?) _ 



[i ?] ^ - C/v ^ [(M*^ -'1 Z 

[r^M) s 9 ] C". t u (tTck) V*- j 






Thus TTo. and ( p J , can be found through known or assumed 

TK. 

quantities • 



Determination of T 3h’ T 5h 

(Tj^ -7~6k.^ 

fak -U) 



E — 

rec 



C c (T;k ~ <2k^) 



where 



Ci. = C\ ^ 

c, = c<v. 



k -f) 



therefore , 



TsL = f [' t (‘ =/A )-<9](u -T(k) 



T 2 ft and can be found once ~iT c ^ ^- s determined. Assume 



5h 



that 



+ 1 00°R 



2. The Rankine Cycle . This is a bottoming steam cycle that 



recovers heat from the closed Braytcn cycle that would other- 
wise be discharged overboard through the cooler. The rela- 
tively hot gas leaving the regenerator is pre-cooled in the 
Waste Heat Boiler (WHB) before being sent to the cooler which 
brings the temperature down to compressor inlet conditions. 

A T-s diagram of the Rankine cycle is shown in Figure 5 * 



T* 




> 

S 

Figure 5 - 
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Major components of the cycle are the WHB, condenser, 
steam turbine, and feed pump. It is assumed that the feed 
pump is driven with power from the turbine. 

The analysis was conducted by assuming values of T^ s , 

AT, ^p( ,5 )(pump isentropic efficiency), (turbine isen- 

tropic efficiency), and 7 WHS (waste heat boiler efficiency). 

The following are defined : 

enthalpy at point is 

absolute temp at point is 

pinch point on Brayton side 
saturation point on steam side 
mass flow ratio (m^/m^) 

quality of steam exiting turbine 
absolute pressure at point is 

power (net) from the steam cycle 

thermal efficiency of the combined cycle 
mass flow of steam in Rankine cycle 

Rankine cycle parameters were chosen for maximum woric . 

The procedure was as follows: 

(1) Compute T^ from the relation 

T 33 = T 6 - AT 

(2) Assume a value for T__. T. < T < 662. 0°F 

PP 1 PP 



T • 
is 

PP 

SF 

z 3 

X 

P is 

*R 

% 

*s 
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(3) Compute T SP from the relation 

t sf = T pp - AT 

(4) Obtain Pg AT * h SF> h^ s from steam tables. 

(5) Compute hji^g from the relation 

h 4 s =h 3 S~ ^ s H h 3 S” h 4 S 60 ) ge* h^ s (/s ) from steam tables 

(6) Obtain x (turbine exit), h^ from steam tables 

(7) Compute from the relation 



- L + dhz 

is ' ?p( ,0 



(8) Compute jS from the relation 




(9) Compute W p from the relation 




*(10) Iterate with T (step 2) until W p is maximum. 



(11) Compute Ty from the relation 



Cf( T,)-T 7 = C?(T n ) -Tfp -k is )/p»6 



(12) Calculate 9-f from the relation 
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III. VOLUME ANALYSIS 



1 . Turbomachinery . An analytical procedure for the minimum 
volume design of turbomachinery was developed by Lee (ref. 
10). The following equations and mathematical expressions 
summarize that portion of the procedure which are applicable 
to compressor and turbine designs in this thesis. Volume 
obtained by Lee's method neglects the volume of casing, inlet 
and outlet plenums, and bearing assemblys. 

a. Turbine . Free vortex blading and axial exit stages 
were assumed throughout the analysis. 



The following are defined t 

S.F. - factor of safety 

U + - max. allowable tip speed 



C L 

V* 

If 

h 

k 



a 
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flow coefficient at tip 

hub-tip ratio at inlet 

hub-tip ratio at exit 

clearance factor 

specific volume at inlet 

max tip radius 

polytropic coefficient 

isentropic exponent 

density of blade material 

stress concentration factor at root 

amplification factor 

correction factor for taper 
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dist. from tip to 1st principle axis 

chord 



-Kzi ~ normalized moment of inertia 

Vol^ - volume of turbine 

The volume was computed by assuming values for SF, U^, 0^, 
f Cj^ i K » j)^ , • a i , Kjj t and 0 • 

r^ t n, and yfi are computed from the following relations: 

- (m e ’</!>, a, .irfi-A 1 )])^ 

= K 

Compute , C 2 » C^# from the following: 

c, = o.z(a*i) 14'' up r/ 

C i = z(A-£1Z.)/(SF •«.) 

c 3 = z&t k (i-y«y"/(sF-<) 

Ca.~ ~ Uif- /z 

where A and B are material dependent constants. The following 
can then be calculated: 

« =C,[( 
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d = c z * c 3 (i u(j-r) 

x>N * +z)f)^ .)))-(', -y;T^*y5 



J>D = ZV(l-y i y''(n-i) C s - Z)>C 4 



The axial chord for a stage is given by b^. For a 
large number of stages, summation of b^ for the whole turbine 
can be approximated by an integral: 








* 



^ . 



where 

J(k) = i (%) A •'O/f] ^ 

This integral can be evaluated numerically thus, 



Vol r = zirr t L c L I T 

b. Compressor . An approach similar to that of the 
turbine was used for the axial compressor. Axial exit stages 
were assumed in the analysis. 

The volume was computed by assuming values for SF, U^., 
C-^, K, jP m> e< , a, x 1 , K I( , Q, and (blade material 
yield stress) . 
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r^, n, and V/3 are computed from the following relations: 

r* r (^s V* /[ tifTT (\ ~ /^)] ) 

K 

n = K -(K~')/9cfo 

^ =[' - - *>‘]* 

Compute C 1 , C2* and from the following: 

C, = X./kj, 

Ci = -Q p fy)0.f / 2 

c, = 2 Ty/fc ' ^ 

Then calculate the following: 

W = C, [('/ ' yj/0 t >'j ! (V r ‘ ".■?'? ^)/(0f 2 r ■ 7«5 
D - Cj ’■C-O-V 2 ) 

dm = c.r-^/r.fty 2 ;/^ 2 '■ 7 2zzf)(i^j)'' 

- (l-1>)0 (fS + /1<?S *.722S»y 

i- .1? y^-y) (i'-v)’ 1 /(p+ 2 + .izis^) 

- (0r Z +.722Sf) z ] 

J>JD = -ZCtV 
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As with the turbine define 



I 



T 




where 



-x>i>-n)/y] c\m 



This integral can be evaluated numerically thus 



voi c = Zirq. l C L l T 



2. Heat Exchangers . The volumes of the Regenerator, Re- 
cuperator, Cooler, and Heater were computed using the method 
of minimum volume design developed by Lee (ref. 10). In 
essence Lee found that the optimal diameter ratio (D^/D) for 
these heat exchangers can be established as a simple relation 

to the pressure ratio (7T C ) of the cycle. The results of the 

AP 

optimally distributed pressure drops ( ) as reported by 

Lee were also studied. Parametric equations relating pressure 
drop to pressure ratio were developed and used in this analysis. 
No consideration has been made on computing the volume of the 
shell, headers, or other appendages. Lee's relations were 
developed from basic heat transfer and pressure drop relation- 
ships . 
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a. Regenerator . The regenerator is assumed to be a 



counter flow shell and tube heat exchanger. High pressure 
is inside the tubes as is the case for current design prac- 
tices . 



The following are defined : 



(—) 

• p 

Bh.) 

/Recj 



Ah 

D 

AT 



i Rwj 



f 

/ ( 

k 



% pressure drop across the regenerator 
(both steams) 

diameter ratio for regenerator 

enthalphy difference across the interior 
stream of Reg. 

tube inside diameter 

log mean temp, difference of hot and cold 
streams 

Prandlt number 
density 

dynamic viscosity 
thermal conductivity 



/AP\ /djA 

Values for \p~J^eg and Re^ can be es ' t ^ nia ' ted from 
the following relations: 



(T^Reg = - - 0284Tr c < f ° r (¥]r = -1. Er = -88) 

(-§) R eg ” -557 * .0933 T c 
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Subscript 1 refers to quantities inside the tubes 
whereas subscript 2 refers to quantities outside the tubes. 
A Re g can be computed as follows: 




where 



_ . , -(.ti i.22 _ -,£64- -.4-1 

T- 17. 2.1 4 k, Mi 'r, 



y - j> u * et ,A{ 



The volume can then be computed as follows: 



Vol Reg ~ A Reg l p /Reg 



A? 



-0.4.1 



b. Recuperator . The recuperator is similar in design 
to the regenerator only it refers to the heater cycle. 



The following are defined : 







( p ) 


Rec 


/-M 




U / 


Rec 


/AP\ 




(TV 


hhr\ 


( 4 - p ) 


PFB 


Ak, 


CO 

o 



% pressure drop across the recuperator 
(both streams) 

diameter ratio for recuperator 

% pressure drop across the heater side 
(outside tubes) 

% pressure drop across the PFB 

Enthalpy difference across the interior 
stream of recup. 
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can be estimated from the 



/AP\ / j > l \ 

Values for(— y Rec and^—^ 

following relations 

m?\ - ^ \ 

(, P ^Pca l P /Tb V r y/iTr^ ( P /PF6 




Subscripts 1 and 2 refer to quantities inside and outside 
the tubes respectively. A Rec is computed as follows: 



a*. •(%) r, v, -v ai, " ( i . &),.,r ] 






where V and Y are as previously defined. 
Then compute the volume 



Vol = \ /-Af \ 

VOi Rec a Rec ( P /Rec 



c. Cooler . The cooler is assumed to be a shell and tube 
heat exchanger with the gas inside the tubes and salt water 
outside . 



The following are defined : 



(f)cir 

(~t/ clr 

Ah, Clr 



v 

V WKB 



2 



% pressure drop across the cooler (air side) 

diameter ratio for cooler 

enthalpy difference across the interior 
stream of cooler 

sp. vol. of working fluid at WHB exit 
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% pressure drop across the WHB (air side) 

% pressure drop across the heater 
(inside tubes) 




polytropic cooling water pump efficiency 



V.. - sp. vol. of cooling water 



(& p ) (2k\ 

Values for \ p / C lr andVjo/Clr are estimated from the 
following relations: 





= -omr, * .076 

Subscripts 1 and 2 refer to quantities inside and outside 
the tubes respectively. A clr coin P u ‘ te<i as follows: 



Ac, 



Akutr 



i.H 




where v , Y are as previously defined and 
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The volume is thus 

,, , _ , /AP\ -.41 

°^clr ~ A clr ( P /clr 



d. Heater . The heater is assumed to be a shell and 
tube heat exchanger with the working gas in the tubes and 
combustion gas out. 



The following are defined : 

- diameter ratio for heater 



) htr 
4h, htr 



enthaphy difference across the interior 
stream of heater 



Values of (4%; ri . (^)htr 2 , aRd <22>> 
estimated from the following relations: 



are 



htr 






6 ? 



\ p J hfr . ~ l P Ati- 



,/AP\ f^\ 

• 1 fi (—Air 



where 



- US 



/* ’(*/«>)" ( a /a) 



2- 




■A) 



(' 



lh)krr = -^ ?7r c + .2/+ 



Subscripts 1 and 2 refer to quantities inside and outside 
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the tubes respectively. Compute A htr as follows: 





r*i 







where , X , °<- are as previously defined. 



Then compute Volj^ r 




e. Waste Heat Boiler. The boiler is assumed to be a 



once through cross-counterflow type, as shown in Figure 6. 
High pressure is inside the tubes thus saturation pressure 
must be greater than the working pressure of the closed 
Brayton cycle . 




.water in 



steam out 



Figure 6. 
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Boiler volume is estimated using heat transfer relations 
and volumes from similar U.S. Navy COGAS (Combined Gas Turbine 
and Steam) designs (ref. 11). 



The following are defined : 



C ps 


- 


ave. specific heat on 


steam side 


C P 


- 


ave. specific heat on 


air side 


NT'U 


- 


Number Transfer Units 




#Ps 


- 


number of passes 




a kt 


- 


Heat transfer area 




U 


- 


overall heat transfer 


coefficient 


V T 


- 


total volume 





From a DTNSRDC (David Taylor Naval Ship Research and 
Development Center) study (ref. 11) the following relations 
are obtained: 

= CMJJl = C^s 0r rnc C P 

CliftX rn s C 

EK =(/»/,£' C?: ' - 0 

Mr U =(-‘‘r , s)/n{i '(/<*){>> [ I - C.<(| - Ee. /tp ‘)/(cf . j 

AfTU * CtilN 
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thus 



For similar designs assume V T 



U 



I 






U 



II 



Vj CMINj NTUj 

VT^ * CMIN^ * NTU^ 




^CMINjj NTU.j ) 



CMIN t NTU t 



The term/V 



II 



can be evaluated from known data 



.CMINjj NTUjj 



on previous designs to give 



V WHB = 11.53 * CMIN * NTU 



But it should be noted that the above equation assumes 

Mp \ 

constant '~J~ > ^3 and boiler power. If as before Vol A.^ 
then from heat transfer 

A ht * f * 

thus 

V °< ty\ s Aln £ 

Correcting the above equation for power gives 
V WHB = .0017 CMIN * NTU * rn, ~ Ah, 

As stated previously Vol °< 
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Correcting again for pressure drop yields 



A? 

? 



1 



-M 



V 'a j HB = * 000 ^5 CM IN * NTU ~m s * 




CMIN in BTU/s - °R 
m s in lbm/s 



4h s in STU/lbm 

f. Condenser . The condenser is assumed to be a shell 
and tube crossflow heat exchanger with salt water inside the 
tubes and steam outside. 

The following are defined . 

D - inside tube diameter 

N - number of tubes 

L - length of tubes 

g - Gravitational acceleration constant 

- pressure drop on the salt water side 
f - friction factor 

fc - density of cooling water 

V - velocity of cooling water 

A c - cross sect, free flow area on water side 
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w 



c 



q 

c 



c 



AT 

c 



X 



h 



fg 



m 



s 



flow rate of cooling water 
heat flow in condenser 

w {xj specif ic heat of cooling water (C ) 
c p 

temp, difference across the cooling water side 

steam quality entering condenser 

enthapy difference between sat. liquid and 
sat. steam 

mass flow of steam entering condenser 



From basic heat transfer relationships 
A HT = ^ DNL 



where 



N = 



4 A c 



S^P C £ 

^ 2 f p<_ Vc 



therefore , 



ZAc 3 AP C 
A N t' ff c y t i 



where 
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thus 



_ 2 ^ AP t h(, x x 



For similar designs assume V T ** A KT 



thus 




~ I for similar designs assuming same JS P 

V 



fas x) x 

[ M S * )lL 



therefore 

v x 





£ 



The quantity in brackets can be evaluated from known 
data on previous designs (ref. 11) to give 

V cond * 2 ?'3 V ft3 

m s in Ibm/s 
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IV. DESIGN OF A 25,000 HP PLANT 



A proposed design of a 25»000 HP coal-fired combined 

cycle was completed. Coal combustion is accomplished using 

a pressurized fluidized bed (PFB). Design parameters were 

eventually chosen after investigating the effect of changing 

pressure ratios, varying (— p —J T in the Brayton cycle, varying 

E^ (regenerator effectiveness), different compressor inlet temp. 

(T^) and steam condensing temp. (T^), and variations in Brayton 

working pressure (WP) and saturation pressure (P m . ). 

inin • 

a. The Fluidized Bed . The pressurized fluidized bed 
operates at a pressure of 3 to 10 atmospheres. Pressure is 
maintained by an exhaust driven compressor. For proper com- 
bustion incoming air must be preheated by the recuperator to 
at least 700°F. The bed material is primarily limestone which 
serves to capture sulfer in the coal and thus minimize SO2 
emissions. A bed temperature of 1?00°F is generally used to 
keep N0 x emissions at an acceptable level. 

The fluid bed itself is about 4 ft. deep. Preheated 
combustion air enters the bed via a perforated bed plate under- 
neath the bed. Pressure drop across the bed is assumed to be 
10 $. 

The PFB is used to heat the working fluid in the closed 
Brayton cycle. The working fluid leaves the regenerator and 
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enters the tubes in the PF3. The working fluid is first 
heated by convection from the gas leaving the bed and is then 
ducted to tubes submerged in the fluid bed for final heating 
to the turbine inlet temp (T^). (Figure 7) 



exhaust 




Figure 7. 
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b. Assumptions . The computer program was run with the 
following input values: 



% (p) 


= 


.88 


'/f (P) 




.91 


T - e r 




.88 


|0 


= 


.80 


MO 




.88 


(^)r 


S 


.10 


? j'JHB 


= 


• 94 


T 1 


= 


540.0 °R 


T 4 


= 


2060.0 °R 


T 1S 




540.0 °R 


At 




+35.0 °R 


wp 




135.0 psi 


p 

SAT 




200.0 psi (want at least 50 osi greater 
than WP) 


Power 




25000 HP 


T ih 


= 


535.0 °R 


T 7h 


- 


760 °R (limited by acid corrosion in the 
turbine) 


9Jf)k 




.85 




= 


.89 


0 


— 


.90 (11$ excess air) 



In addition it was assumed that scoop injection could be 
used for both the cooler and condenser at speeds greater than 
10 knots, therefore power requirements for salt water circulating 
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pumps were excluded in the calculations. 



c. Results . The computer output is provided in Appendix 
3. Note that the saturation pressure for pressure ratios of 
2 and 3 were below the 200 psia specification. This is due 
to the fact that the low pressure stream leaving the regener- 
ator is at a temperature below the saturation temp, at 200psia. 

vs. il ^ and Vol +rt+r>1 vs. ir rt is shown in Figure 



A plot of ty-f- vs. TT c and Vol^ ^ ^ vs • IT c 



8 . 

Note that the volume increases sharply above TP = 4 
due largely to the rapid growth in the waste heat boiler. 
Maximum occurs at around li c = 6 . To investigate the 

effect of on and Vol^, the program was run for 

-4^) ip = .04, . 06 , .08, .10, and .12 between IT c of 4 and 9* 



The results plotted on Figure 9 show that significant 
gains in J/ r can be realized with very little increase in V T 

for the range of 4-6 pressure ratio. The effect of changing 

/AP \ 

E P (regenerator effectiveness) in this range with a '~p ~ = 
.05 was further investigated. 



The results listed in Table 1 show no increase in effi- 
ciency as is reduced. Furthermore, volume increases sharply 
as a greater percentage of the work is being done by the Rankine 
cycle as E^ decreases. 

Thus for greatest efficiency with minimum volume, cycle operation 
should be at the highest E p practicable. 
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TT C 



Figure 8, 



43 



Figure 9. 
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Table 1 



E r t Vol t (« 3 ) 





.88 


.452 


335 


3 


.80 


.423 


646 




.50 


.402 


7038 




.88 


.464 


4o8 


4 


.80 


.454 


811 




.50 


.432 


4350 




.88 


.474 


724 


5 


.80 


.464 


1083 




.50 


.447 


3741 




.88 


.475 


1050 


6 


.80 


.471 


1396 




• 50 


.456 


3508 




.88 


.473 


1450 


7 


.80 


.471 


1750 




.50 


.460 


34 35 
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To investigate the effect of changes in T^ (compressor 
inlet temp.) runs were made on 4 different temperatures and 
the results shown in Figure 10. As expected efficiency de- 
creased with increasing with very little change in the range 
of interest 4-4.5 TT C * Volume on the other hand decreased 
slightly for TT Q less than 3-5 hut increased significantly 
above that. Thus for best efficiency and least volume T^ 
should be kept as close to cooling water temperature as possible. 

Finally, changes in working pressure of the Bray ton cycle 
were investigated. Decreases in working pressure may require 
additional volume of Brayton cycle components yet at the same 
time allow a decrease in Rankine cycle steam pressure (recall 
50°R pressure differential between boiler streams). This will 
in turn increase cycle efficiency yet further increase volume. 
These results are born out in Figure 11. 

d. System Selection . In accordance with the design 
criteria the highest thermal efficiency considered is 0.48? 
(Figure 11). Thus it is desired to select the plant with min- 
imum volume at = .482. 

Figure 11 shows that only those alternatives with Brayton 
cycle working pressures of ?0psi and lOOpsi meet the efficiency 
requirement. Of these two the lOOpsi case has minimum volume 
at desired efficiency. 

However, the problem does not end here. If one looks at 
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Figure 10. 
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1400 



1000 



Vol t 

600 



200 




3 4 5 -rr 6 7 8 
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the steam quality leaving the turbine at the TT c = 5 » WP = 

100 psi point it is of the order of 0.84. This is a little 
lower than normally desired for design. In fact, the steam 
quality of all alternatives thus far considered within the 
range of interest is low. If a steam quality of order .8? is 
desired the current design pt. can change in one of 3 ways 
as follows s 

1) Raise steam condensing temp (T^ s ) 

2) Lower minimum saturation pressure (? m ^ n ) 

3) or a combination of the above two 

Through trial and error it was found that case 1 gave the 
greatest benefit with the least penalty. Computer runs were 
made in the 4-5 Tf range with the steam condensing temp 
raised to 560°R for both the 70 and lOOpsi alternatives. The 
results show a substantial increase in quality, slight decrease 
in efficiency and very little change in total volume (Table 2). 



(psi) 

Working Pressure 


H c 


n 


X 


Vol t (ft 


3 ) 


70 


4.5 


.480 


.864 


948.3 




70 


5.0 


.482 


.871 


1140.6 




100 


5.0 


.479 


.862 


956.5 






Table 2. 








The 70 psi 


designs have better 


steam quality and < 


h than 


the best 100 psi 


case. In 


addition, 


at 4.5 


IT the volume is 
c 


lower. Between 


the two 70 


psi cases 


above a 


trade-off 


exists 



between y ^ and Vol+. Assume that volume comparable to the case 
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of TT = 5 t WP = 100 psi, = 540°R discussed above (Figure 

11) is desired. If this is the case then the new design point 
is at a Brayton cycle working pressure of 70 psi, JT c = 4.5, 
T^ s = 560°R. The results of this design are summarized in the 
following computer printout: 
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BRAYTON CYCLE: PRESSURE RATIO =4.50 

WORKING PRESSURE = 70.0 



TOTAL PRESSURE 


DROP = .05 






MASS FLOW (LBM/S) =206. 9 3 






STATE TEMP (DEG R) 


SUMMARY ‘.THERMAL EFF=0.480 


1 


540.0 




HEATER EFF=0 .933 


2 


876.0 




SFC ( LBM/HP -HR ) =0 . ■ 


3 


1403-9 




VOL TOTAL ( CU FT ) = 


4 


2060.0 




POWER (HP) =25000.0 


5 


1475-9 




WORK FRAC BY STM 


6 


948.0 






PP 


836.2 






7 


809.2 






COMPONENT 


DPP DHD 


VOL(CU FT) 


, FSGEN 


.010 .950 


26.5 


EFF=0 . 88 


TURBINE 




5.5 


EFF( P=0 .91 


COOLER 


.009 .490 


193-1 




COMPRESSOR 




29.7 


EFF( P) =0 . 88 


HEATER 


.0006 






BOILER 


.030 




EFF=0 .94 



(PHI=0.90) 



HEATER CYCLE: PRESSURE -RAT 10=4 . 79 

TOTAL PRESSURE DROP =.43 
AIR FLOW ( LBM/S ) = 36.85 
FUEL FLOW ( LBM/S )= 3*06 



STATE TEMP (DEG R) 



1H 


535.0 






2H 


835.2 






3H 


1440.5 






4H 


4809.8 






5H 


1503.9 






6H 


945.0 






7H 


760.0 






COMPONENT 


DPP 


DHD 


VOL(CU FT) 


RECUP 


.313 


.650 


1.9 EFF=0 . 91 


TURBINE 






0.7 EFF(P)=0 .89 


HEATER 


.019 


.430 


31 .4 



COMPRESSOR 10.8 EFF(P)=0.85 



RANKINE CYCLE: SATURATION PRESSURE = 120.0 

STEAM FLOW ( LBM/S )= 5.69 
FLOW RATIO =0.0275 



STATE TEMP (DEG R) 


PRESSURE (PS I) 


IS 


560.0 


0.9 


2S 


575.0 


120.0 


SF 


801.2 


120.0 


3 s 


913.0 


120.0 


4s 


560.0 


0.9 QUAL=0 . 864 


COMPONENT 


VOL( CU FT) 


BOILER 


507.9 




TURBINE 


6.5 


EFF( IS ) = 0 . 88 


CONDENSER 


134.2 
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V. APPLICATION 



The application of this system to Navy ships has many 
implications. One of which is the fact that neither Closed 
Bray ton nor combined cycles have been previously used. However, 
there is a considerable effort at this time to get COGAS 
(Combined Gas Turbine and ’Steam) systems on board Navy ships. 
These systems are all limited to burning oil fuels. The Closed 
Brayton combined cycle when equipped with a multi-fuels heater 
offers the additional flexibility of burning the cheapest or 
most readily available fuel. 

Advantages of the Closed Brayton combined cycle ares 

1) The Brayton cycle operates at temperatures substan- 
tially below the material limit than in the open cycle gas 
turbines resulting in cheaper, more reliable, and longer life 
components. 

2) Low pressure steam systems are well proven, and have 
significantly less water chemistry problems than higher pressure 
ones . 

3) As load decreases and fuel flow is cut back, the con- 
trol valves (Figure 12) between the compressor outlet and cooler 
inlet can be opened as necessary (keeping max. temp and speed 
constant) so as to ensure that the vector diagrams of all sec- 
tions of the turbomachinery blading remain unchanged. Turbine 
and compressor efficiencies will remain essentially the same 
producing high 0 ^ at part load. This requires use of a CRP 
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Functional Component Arrangement 



Figure 12 




s-. 

• H 
< 
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Exhaust 



(controllable reversible pitch) propeller to maintain constant 
shaft speed at part loads. 

Disadvantages of the Closed Brayton combined cycle are: 

1) Higher overall system weight and volume due to the 
additional heat exchangers and piping required for the Closed 
Brayton cycle. 

2) A complex, integrated control system is required to 
coordinate control of the combined cycle, fuel flow, and CRP 
at part loads. 

3) Fluidized bed combustion is not a fully developed 
engineering concept. Land based prototypes have been built 
and operated with limited but increasing levels of success. A 
shipboard prototype is perhaps 5 years in the future with pro- 
duction of a narinized version at least 10 years away. 

4) If coal is to be used as the primary fuel, it should 
be noted that it has a lower heating value per ton than oil. 

The result is a larger potentially more expensive (greater 
acquisition cost) ship for the same military mission. 

3) Reaction time with coal fired fluidized bed combustors 
to load changes are generally slower than that for oil heaters 
necessitating in a separate fuel oil system for rapid increases 
in load. 



34 



VI. CONCLUSIONS AND RECOMMENDATIONS 



The Closed Brayton combined cycle offers an attractive 
alternative system for many shipboard applications. As with 
the first Arab oil embargo, subsequent cut-offs of Arab oil 
could result in decreased steaming days and a lower overall 
readiness of the U.S. Navy’s defense posture. 

If our ships could have the option of burning the most 
readily available fuel, this reduction in readiness can be 
avoided by switching to coal, a fuel source in great abundance 
in the U.S. 

Paramount in the engineering of a successful multi-fuels 
propulsion plant is the development of a marinized fluid bed 
combustor. The design of heat exchangers and other turbo- 
machinery used in the combined cycle are well within the state 
of the art of present technology. 

It is thus that the author concludes the following: 

1) The design of a modern coal burning combined cycle is 
within the reach of present technology. 

2) This combined cycle would offer high efficiency at 
both rated load and part load due to the unique controlling 
aspect of the Closed Brayton cycle. 

3) Weight and volume of a Closed Brayton combined cycle 
would be greater than that of an equivalently rated COGAS type 
system. 

4) Coal offers the benefit of a readily available, 
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relatively cheap fuel at the expense of increased ship size 
and acquisition cost due to lower BTU heat value per ton when 
compared to oil. 

The author makes the following recommendations: 

1) Development of a land based test site of a Closed 
Brayton combined cycle equipped with a fluid bed heater should 
be considered. 

2) A pressurized fluidized bed combustor suitable to 
marine applications should be developed. 

3) A shipboard prototype system should be integrated 
on a U.S. Navy ship to attempt to quantitatively define the 
shipboard impact of the system. 
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APPENDIX A 



The computer program performs a thermodynamic analysis 
and component volume computation for various specified 
Brayton cycle compressor pressure ratios. 

The program accepts as input the following quantities: 



7<-(p> 

?+(p> 

Er, 






iS) 



AP 

F/ T 



/ WHB 

m 

x i 



ns 

AT 



polytropic compressor efficiency (Brayton 
cycle ) 

poly tropic turbine efficiency (Brayton 
cycle) 

regenerator effectiveness 

isentropic pump efficiency (Rankine 
cycle ) 

total pressure drop around Brayton cycle 
loop 

efficiency of the waste heat boiler 
compressor inlet temp. (Brayton cycle) 

turbine inlet temp (Brayton cycle) 

steam condensing temp. 

- min. temp diff . between streams in WHB 
(at P? and Tg) 



?S^t - min. saturation pressure in the Rankine cycle 



Power 




7Vf>h 

A* ( is) 



net power out of the combined cycle 
compressor inlet temp. (heater cycle) 

turbine exit temp, (heater cycle) 

polytropic compressor efficiency (heater 
cycle) 

isentropic turbine efficiency (Rankine 
cycle) 
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polytropic turbine efficiency (heater cycle) 
equivalence ratio 

working pressure in the Brayton cycle. 



%(P) K 

0 

WP 

Output is as shown in Appendix B. 

The main program performs the thermodynamic analysis of 
the closed Brayton cycle and the Rankine cycle. Functions per- 
formed by the various subfunctions are as follows: 

1) Subroutine HEATER performs a thermodynamic analysis of 
the heater cycle. 

2) Subroutine DHD computes-^- for the regenerator, re- 
cuperator, cooler and heater. 

. hP 

3) Subroutine DPP calculates ~y ~ for the regenerator, re- 
cuperator, cooler, heater and boiler. 

4) Subroutine TURVOL finds turbine volumes. 

5) Subroutine COMPVOL finds compressor volumes. 

6) Subroutine REGEN computes the regenerator volume. 

7) Subroutine RECUP computes the recuperator volume. 

3) Subroutine COOLER computes the cooler volume. 

9) Subroutine COMBUSTOR computes the heater volume. 

10) Subroutine WHBVOL computes the waste heat boiler volume. 

11) Subroutine CONDVOL computes the condenser volume. 

12) Subroutine STMTAB is a computerized version of the steam 
tables. Reproduced from ref. 15* 
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DIMENSION CR(30) ,SATLIQ(3) ,STEAM(3) ,TPP(12) ,WR(1 1) 

REAL NCP,NTP,NPIS,NTIS,NWHB,NT,Ll,L2,KO 
REAL NCPH,NTPH,MB ,MS,MA,MF ,NH,N , IT,K II ,NPP 
REAL LHV,M,M1,M2,NTU 

COMMON /STM/PRESS , TEMP ,QUAL , SATL IQ, STEAM 
COMMON /HTR/NH ,PRH ,DPPH ,T2H , T3H , T4H , T5H , T6H ,MA ,MF 
COMMON /DIA/DHDREG,DHDCLR,DHDHTR,DHDREC 
COMMON /PRES/DP PWHB ,DPPREG,DPPCLR,DPPREC ,DPH1 ,DPH2 ,DHDHT 
COMMON /COMB/DPPRG,DPPRC,DPPHE,DHDHR,DPHE2 
READ (5, 10) NCP ,NTP ,ER,NPI S ,NTI S ,DP,NWHB ,T1 ,T4 ,T1S ,DT ,PMI N,K 
10 FORMAT (7F4.3,5F7.1,I2) 

READ (5, 15) PWR,T1H,T7H,NCPH,NTPH,PHI,WP 
15 FORMAT (F8.1,2F6.1,3F5.3,F6.1 ) 

READ (5 ,20) ( CR ( I ) ,1=1, K) 

20 FORMAT (<K>F6.2) 

R= .068544 
DO 200 1=1, K 
PR=CR( I ) 

Z=CP1(T1)+(R/NCP)*L0G(PR) 

T2=EXP(Z/.2475) 

CALL ITER (T2,Z) 

TR=1/(1-DP)/PR 

Z=CP1(T4)+R*NTP*L0G(TR) 

T5=EXP(Z/.2475) 

CALL ITER (T5,Z) 

IF ( T2 .GT. T5) GO TO 30 
T6=T5-ER*(T5-T2) 

T3=T2+ER*(T5-T2) 

GO TO 40 
30 T6=T5 

T3=T2 
40 L1=T1 

L2=1122.0 

IF ( L2 .GT. T6) GO TO 50 
GO TO 60 

50 L2=T6-( T6-T1) /10. 

60 TPP( 1 )=L2 

WR1=0. 

C=(L2-L1)/10. 

DO 70 J=l,ll 
TSF=TPP( J )-DT 
T3S=T6-DT 
T=TSF-460 

CALL STMTAB (2,T,0.) 

PSF=PRESS 

HSF=SATLIQ(2) 

T=T3S-460 

CALL STMTAB (7,PSF,T) 

H3S=STEAM(2 ) 

S3S=STEAM(3) 

T=TlS-460 



61 



CALL STMTAB (2,T,0.) 

P4S=PRESS 
V 1S=SATL I Q ( 1 ) 

H1S=SATLIQ(2) 

S1S=SATLIQ(3) 

CALL STMTAB (9,P4S,S3S) 

H4SIS=STEAM(2) 

H4S=H3S-NTIS*(H3S-H4S IS) 

CALL STMTAB (8,P4S,H4S) 

X=QUAL 

V4S=STEAM(1) 

H2S IS=H1S+V 1S*(PSF-P4S ) 

H2S=H1S+(H2SI S-H1S) /NP I S 
T2S=T1S+15. 

B=NWHB*(CP2(T6)-CP2(TPP (J) ) )/ (H3S-HSF ) 

Z=ACP (TPP(U ) )*TPP(J )-B*(HSF-H2S ) /MWHB 
T7=TPP(J ) 

CALL ITER2 (T7,Z) 

WR(J )=B*( H3S+H1S-H4S-H2S) 

IF ( J .EQ. 1) GO TO 68 
IF ( T7 .LT. Tl) GO TO 75 ' 

IF(PSF .LT. PM IN) GO TO 75 
IF ( WR ( J ) .GT. WR1) GO TO 68 
I F (X .GT. .87) GO TO 75 

68 WR1=WR( J ) 

TPP(J+1) =TPP(J )-C 

70 CONTINUE 

GO TO 76 

75 IF(ABS(TPP(J )-TPP(J— 1) ) .LT. .1) GO TO 76 

L1=TPP( J ) 

L2=TP P ( J - 1 ) 

GO TO 60 

76 WB=CP2(T4)+CP2(T1)-CP2(T5)-CP2(T2) 

QIN=CP2( T4)-CP2(T3) 

NT=(WB+WR(J ) )/QIM 
MB=.707*PWR/(WB+WR(J) ) 

MS=B*MB 

QB=MB*QIN 

CALL HEATER (QB ,T1H ,T7H ,NCPH,NTPH,PHI,T3 ,T4) 

CALL DHD(PR) 

DHDHT=DHDHTR 

K=1 

77 CALL DPP(T6 ,T7 ,PR,DP,ER,T2,T3 ,DPPH,T4H ,T5H ,T4 ,PRH,WP ,MB ,MA 
* ,MF ,K) 

G=1.33 

N=G/ ( G— ( G— 1 . )*NTP ) 

P=WP* (PR-DPPREG) 

V=1./ARH0(T4,P) 

V MNTURB=0.0 

CALL TURBVOL (N,MB ,V,TR,T4 ,VMNTURB ) 

G=1 .39 
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N=G/(6-(G-l.)*NTPH) 

P=14. 7* (PRH-DPPREC) 

V=1./ARH0(T6H,P) 

M=MA+MF 

VHTRTURB=0 .0 

TRH=1 . / (1 . — DPPH) /PRH 

CALL TURBVOL ( N,M,V ,TRH,T6H ,VHTRTURB) 

T=T1S/T3S 

P=P4S/PSF 

N=1 . / (1 . -LOG(T) /LOG(P ) ) 

VSTMTURB=0.0 

CALL TURBVOL (N,MS,V4S ,P ,T3S ,VSTMTURB) 

6=1.4 

N=G/(G-(G-1. )/NCP) 

V=1./ARH0(T1,WP) 

VMNC0MP=0.0 

CALL COMP VOL (N,MB,V,PR,VMNCOMP) 

V=1 . /ARHO (535. ,14.7) 

VHTRC0MP=0.0 

CALL COMPVOL (M,MA,V ,PRH,VHTRCOMP ) 

D= .125 
VREGEN=0.0 

CALL REGEN (T2,T3 ,T5 ,T6 ,WP ,D ,MB ,PR ,DPPREG,DHDREG, VREGEN) 
VRECUP=0.0 

CALL RECUP (T2H,T3H,T5H,T6H, PRH, DPPH, D,MA,MF,DPPREC,VRECUP) 

D= .25 
VCLR=0.0 

CALL COOLER (T7 ,T1 , WP ,D ,NTP, MB ,DHDCLR,DPPCLR ,VCLR) 

DPPRG=DPPREG 

DPPRC=DPPREC 

DPPHE=DPPH 

DHDHR=DHDHTR 

DPHE2=DPH2 

VHTR=0.0 

CALL COMBUSTOR (T3,T4,T4H,T5H,WP,PR,D,MB,PRH,MA,MF,VHTR) 
VWHB=0 .0 

CALL WHBVOL (MB ,MS,T6 ,T7 ,NWHB ,H2S ,H3S ,DPPWHB ,VWHB) 

VC0ND=0 .0 

CALL CONDVOL (MS,X,VCOND) 

VOLT=VMNTURB+VHTRTURB+VSTMTURB+VMNCOMP+VHTRCOMP+VREGEN+VRECUP 
* +VCLR+VHTR+VWHB+VCOND 
IF(K .EQ. 1) GO TO 78 
IF (VOLT .GT. V0L1) GO TO 81 
I F ( (DPPCLR- .002) .LT. 0.0) GO TO 81 
78 K=K+1 

V0L1=V0LT 
V1=VCLR 
V2=VWHB 
GO TO 77 
81 V0LT=V0L1 

VCLR=V 1 
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VWHB=V2 

WRITE(6 ,82) PR, WP, DP, MB 

82 FORMAT( ' 1 ' ,'BRAYTON CYCLE: PRESSURE RATIO ' ,F5 .2/1X, 'WORKING 

* PRESSURE 2 ' ,F5 .1/1X, 'TOTAL PRESSURE DROP 2 ' ,F3.2/1X, 'MASS FLOW 

* (Ibm/s) 2 ' ,F6 .2//10X, 'STATE TEMP (deg R) 1 ) 

WRITE(6 ,83) T1 

83 FORMAT (12X, ' 1 ' ,6X,F6.1) 

WRITE(6 ,84) T2 

84 FORMAT ( 12X , ' 2 1 ,6X,F6.1) 

WRITE ( 6 ,85) T3 

85 FORMAT (12X, '3' ,6X,F6.1) 

WRITE(6 ,86) T4 

86 FORMAT (12X , '4' ,6X,F6 .1 ) 

WRITE(6 ,87) T5 

87 FORMAT (12X , ' 5' ,6X,F6 . 1 ) 

WRITE ( 6 ,88) T6 

88 FORMAT (12X, ' 6 1 ,6X,F6 .1 ) 

WRITE(6 ,89) TPP(J) 

89 FORMAT (12X, ’ PP 1 ,5X,F6 .1) 

WR I TE ( 6 ,90) T7 

90 F0RMAT(12X, ’7' ,6X,F6 .1//5X, 'COMPONENT DPP DHD VOL(cu ft) ' ) 
WRITE (6 ,91) DPPREG,DHDREG,VREGEN,ER 

91 F0RMAT(7X, 'REGEN' ,4X,F4.3,1X,F4.3,1X,F7.1,2X, ' EFF= ' ,F4.2) 
WRITE (6, 92) VMNTURB ,NTP 

92 F0RMAT(6X, 'TURBINE' ,15X,F5.1,2X, 'EFF(P)='F4.2) 

WRITE (6 ,93) DPPCLR,DHDCLR ,VCLR 

93 F0RMAT(6X, 'COOLER' ,4X, F4 .3 ,1X, F4 .3 ,3X, F5 .1 ) 

WRITE (6,94) VMNCOMP ,NCP 

94 F0RMAT(4X, 'COMPRESSOR' ,14X,F5.1,2X, 'EFF(P) 2 ' ,F4.2) 

WRITE (6,95) DPH1 

95 F0RMAT(6X, 'HEATER' ,4X,F5. 4) 

WRITE (6,96 ) DPPWHB ,NWHB 

96 F0RMAT(6X, 'BOILER' ,4X,F4.3,15X, ' EFF 2 ' ,F4 .2/ ) 

WRITE (6 ,100) PRH,DPPH,MA,MF 

100 FORMAT( IX, 'HEATER CYCLE: PRESSURE RATIO 2 ' ,F5 .2/1X, 'TOTAL 

* PRESSURE DROP 2 ' ,F3 .2/1X , 'AIR FLOWClbm/s) 2 ' ,F6.2/1X, 'FUEL 

* FLOW (lbm/s) 2 ’ ,F6.2//10X, 'STATE TEMP (deg R) ' ) 

WRITE(6 ,101) T1H 

101 F0RMAT(12X, ' 1 h ’ ,5X,F6.1) 

WRITE(6 ,102) T2H 

102 F0RMAT(12X, ' 2 h ' ,5X,F6.1) 

WRITE { 6 ,103) T3H 

103 F0RMAT(12X, '3h' ,5X,F6.1) 

WRITE (6 ,104) T4H 

104 F0RMAT(12X, *4h ' ,5X,F6.1) 

WRITE ( 6 ,105) T5H 

105 F0RMAT(12X, ' 5 h ' ,5X,F6.1) 

WRITE (6 ,106) T6H 

106 F0RMAT(12X, '6h' ,5X,F6.1) 

WRITE (6 ,107) T7H 

107 F0RMAT(12X, '7h' ,5X, F6 .1//5X, 'COMPONENT DPP DHD VOL(cu ft) ' 
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* ) 

ERH= (T3H-T2H ) / (T5H-T2H ) 

WRITE (6, 110) DPPREC ,DHDREC ,VRECUP,ERH 

110 FORMAT! 7 X, 'RECUP ' ,4X,F4.3,1X,F4.3,1X,F7.1,2X, ' EFF= ' ,F4.2) 
WRITE (6 ,111) VHTRTURB ,NTPH 

111 FORMAT! 6 X, 'TURBINE' ,15X, F5 .1 ,2X, 'EFF(P)= ' ,F4.2) 

WRITE (6 ,112) DPH2,DHDHTR,VHTR 

112 F0RMAT(6X, 'HEATER' ,4X,F4.3,1X,F4.3,3X,F5.1) 

WRITE (6 ,113) VHTRCOMP ,NCPH 

113 F0RMAT(4X, 'COMPRESSOR' ,14X, F5 .1 ,2X, 'EFF(P ) = ' ,F4 . 2/ ) 

WRITE (6 ,120) PSF ,MS ,B 

120 F0RMAT(1X, 'RANKINE CYCLE: SATURATION PRESSURE= ' ,F6 .1/1X, 

* 'STEAM FLOW(lbm/s)=' ,F6.2/1X, 'FLOW RATIO= ' ,F6 .4//10X, 'STATE 

* TEW (deg R) PRESSURE(psi) ' ) 

WRITE(6 ,121) T1S ,P4S 

121 F0RMATU2X, 'Is' ,5X,F6.1,8X,F6.1) 

WRITE (6 ,122) T2S,PSF 

122 F0RMAT(12X, '2s' ,5X,F6 .1 ,8X,F6 .1) 

WRITE(6 ,123) TSF ,PSF 

123 FORMAT! 1 2X, ' SF ' ,5X,F6.1,8X,F6.1) 

WRITE (6 ,124) T3S ,PSF 

124 FORMAT! 12X, '3s' ,5X,F6.1,8X,F6 .1) 

WRITE (6 ,125) T1S ,P4S ,X 

125 F0RMAT(1 2X, '4s' ,5X,F6.1,8X,F6.1,2X, ' Q UAL= ' ,F5 .3/ ) 

WRITE (6 ,128) VWHB 

128 FORMAT! 5X, 'COMPONENT VOL(cu ft) '/6X, 'BOILER' ,5X,F6 .1) 

WRITE (6, 129) VSTMTURB ,NTI S 

129 F0RMAT(6X, 'TURBINE' ,5X, F5 .1 ,2X, ' EFF( IS)= ' ,F4 .2) 

WRITE(6,130) VCOND 

130 F0RMAT(5X, 'CONDENSER' ,4X,F5.1/) 

WRITE (6 ,135) NT 

135 FORMAT! '1' .'SUMMARY: THERMAL EFF=',F5.3) 

WRITE(6,136) MH.PHI 

136 F0RMAT(9X, 'HEATER EFF=' ,F5. 3, ' (PHI= ' ,F4 .2, ' ) ' ) 
SFC=MF*3600./PWR 

WRITE(6 ,137) SFC ,VOLT,PWR 

137 FORMAT (9X, 'SFC(lbm/HP-hr)=' .F5.3/9X, 'VOL TOTAL(cu f t)= ' ,F6 . 1/ 
*9X, 'POWER(HP)*' ,F7.1) 

RANKIME=WR(J)/WB 
WRI TE( 6 ,138) RANKINE 

138 FORMAT (9X, 'WORK FRAC BY STM CYCLE= ' ,F5 .3 ) 

200 CONTINUE 

. STOP 
END 

FUNCTION CP1(T) 

CP1=.2475*L0G(T)-(3.759E-5*T)+(2.553E-8*T**2)-(4.103E-12*T**3) 

RETURN 

END 

FUNCTION CP 2 ( T ) 

CP2=.2475*T-(1.879E-5*T**2)+(1.702E-8*T**3)-(3.078E-12*T**4) 

RETURN 
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END 

FUNCTION AMU(T) 

TK=5.*T /9 . 

AMU=TK**1 .*5*l.E-6/ (TK+110.4) 

RETURN 

END 

FUNCTION ARHO( T,P) 

R=53.35 

ARH0=P*144 ./R/T 

RETURN 

END 

FUNCTION AK(T,P) 

A= .6325E-5 

B=245.4 

C=12. 

TK=5.*T /9 . 

AK0=(A*TK**.5)/(1.+B*10.**(-1.*C/TK)/TK)*418.6 

AK=( AKO+( 1.464E-5* (16. 02*ARHQ(T,P))**1. 23))*. 578/3600. 

RETURN 

END 

FUNCTION APR(T,P ) 

CP=ACP(T) 

XMU-AMU(T) 

XK=AK( T,P ) 

APR=CP*XMU/XK 

RETURN 

END 

FUNCTION A$IGMA( T,P) 

ASIGMA=17.214*(6228.*AK(T,P))**(-1.41)*(1.46*AMU(T))**1.23* 

*(1./APR(T,P)**.564)*(1./(16.02*ARH0(T,P))**.41) 

RETURN 

END 

FUNCTION ACP(T) 

ACP=. 2475-3. 759E-5*T+5.106E-8*T**2-1.231E-11*T**3 

RETURN 

END 

SUBROUTINE ITER (T,X) 

10 T1=T 

T=EXP( (X+.2475*L0G(T)-CP1( T) ) /. 247 5) 

IF ( ABS ( T -T1 ) .LT. l.E-2) GO TO 20 
GO TO 10 
20 RETURN 
END 

SUBROUTINE ITER2(T,X) 

10 T1=T 

T=X/ACP(T) 

IF ( ABS ( T-Tl ) .LT. l.E-2) GO TO 20 
GO TO 10 
20 RETURN 
END 
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SUBROUTINE HEATER (QB ,T1H,T7H,NCPH,NTPH,PHI ,T3,T4) 

REAL LHV,NH,MF,MA,NCPH,NTPH,N 

COMMON /HTR/NH,PRH,DPPH,T2H,T3H,T4H,T5H,T6H,MA,MF 

LHV=12900. 

FA= 0922 

H7H=ACP(T7H )*T7H 
H1H=ACP(T1H)*T1H 

NH=l-( (1+FA*PHI )*H7H-H1H ) /(FA*PHI*LHV) 

MF=QB/(NH*LHV) 

MA=MF/(FA*PHI) 

T2H=0 .0 
T6H=0.0 

10 CPT=ACP( (T6H+T7H)/2. ) 

CPC=ACP( (T1H+T2H )/2 . ) 

N= ( 1+FA*PH I ) *CPT*NTPH* ( CPT*T7H-CPC*T1H ) 

D=CPC*T1H* (CPC/NCP H-CPT*NTPH) 

PRH=(N/D)**(NCPH/CPC) 

TH2=T1H*PRH**(CPC/NCPH) 

N= ( 1+FA*PH I ) *CPT*T7H+CPC*T1H* ( PRH** ( CPC/NCPH )- 1 ) 

D= ( 1 +F A*P H I ) *CP T*T7 H*PRH**(CPT/NTPH) 

DPPH=1-(N/D)** (1 ./( CPT*NTPH) ) 

TRH=(1-DPPH)*PRH 
TH6=T7H*(TRH)**(CPT*NTPH ) 

I F ( ABS ( T2H— TH2 ) .LT. .1) GO TO 20 
IF ( ABS ( T6H— TH6 ) .LT. .1) GO TO 20 
T2H=TH2 
T6H=TH6 
GO TO 10 
20 T2H=TH2 
T6H=TH6 
T5H=T3+100. 

T4H=T5H+QB/(ACP (T5H/2+1300. )*(MA+MF) ) 
T3H=T2H+(1+MF/MA)*(T5H-T6H) 

RETURN 

END 

SUBROUTINE DHD(PR) 

COMMON /DIA/DHDREG,DHDCLR,DHDHTR,DHDREC 

DHDREG=.557+.0933*PR 

IF (DHDREG .GT. .95) DHDREG= .95 

DHDCLR=.076+.092*PR 

DHDHTR= .214+ .048*PR 

DHDREC= .65 

RETURN 

END 

SUBROUTINE DPP(T6,T7,PR,DP,ER ,T2,T3,DPPHJH4,TH5,T4,PRH,P,M,MA 
*,MF,K) 

REAL M,MA,MF,M1,M2 

COMMON /PRES/DPPWHB ,DPPREG,DPPCLR ,DPPREC ,DPH1 ,DPH2 ,DHDHT 
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DTG=T6 -T7 

DPPWHB=FL0AT(K)*.002 
DPPREG=( ,135-.0284*PR)*DP/.l 
IF (DPPREG .LT. .01) DPPREG=.01 
IF(T2 .EQ. T3 ) DPPREG=0.0 
DPH2=( .0108*PR-.01)*DP/.l 
T1=(T3+T4)*.5 
TT=(TH4+TH5)* .5 
P1=P*(PR-DP) 

P2=14.7*(PRH-DPPH) 

M1=M 

M2-MA+MF 

XMUl-AMU(Tl) 

XMU2=AMU(TT) 

XRH01=ARH0(T1,P1) 

XRH02=ARH0(TT P2) 

B»(M1/M2)**(-1.75)*(XMU1/XMU2)**(-.25)*(XRH01/XRH02)*(P1/P2) 

DPH1=DPH2*DHDHT**3/B 

DPPCLR=DP-DPPWHB-DPPREG-DPH1 

DPPREC=DPPH-DPH2- .1 

RETURN 

END 

SUBROUTINE TURBVOL (G,M,V ,TR,T,VOL) 

REAL M,N,KI1 ,IT 

SF=4.0 

UT=850. 

PHIT=0.6 
VA=0.9 
CL=2.0 
P 1=3 .14159 
RH0M=15.12 
A=9 .63E6 
B=-95. 14 
ALPHA-1.2 
Al-2.5 
Xl= .307 
KIl=3.39E-3 
THETA=0.8 

RT= SQR T ( M*V / ( P H I T*U T*P I' * ( 1 . - V A**2 ) ) ) 
VB=SQRT(1.-TR**(1./G)*(1.-VA**2)) 

Cl = 0 . 2* ( A1 + 1 . ) *1 . /V* ( 1 . - V A**2 ) *U T**2*RT**2 *X 1 /K 1 1 
C2=2.*(A-B*T)/(SF*ALPHA) 

C3=2 . *B*T*(1 . -VA**2 ) ** (G- 1 . ) / (SF*ALPHA) 

C4=-l .*THETA*RH0M*UT**2/2. 

H=(VA-VB)/20. 

SUM-0.0 

X1=0.0 

Y1=VA 

DO 100 J-1,21 

N— Cl* ( (1.-Y1)/(1.+Y1))*PHIT**2+Y1**2 
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D=C2+C3* ( 1 . -Y 1**2) ** ( 1 . -G )+C4* ( 1 . -Y 1**2 ) 
DN=C1*(2.*Y1*(1.-Y1)/(1.+Y1)-(PHIT**2+Y1**2)/(1.+Y1)-(1.+Y1) 
*** (-2)* ( 1 . -Yl)* (PHIT**2+Y 1**2) ) 

DD=2 . *Y 1* ( 1 . -Y 1* *2 ) ** ( - 1 . *G ) * ( G- 1 . )*C3-2.*Y1*C4 
X 2= . 5* ( 1 . / ( N /D ) ** . 5 ) * ( (DN*D— DD*N ) /D**2 ) 

I F ( X 1 .EQ. 0.0) GO TO 50 
SUM=SUM+ (X1+X2) /2 . 

50 X1=X2 

Y1=Y1-H 
100 CONTINUE 
IT=H*$UM 

V0L=2.*PI*RT**2*CL*IT 

RETURN 

END 

SUBROUTINE COMPVOL (G,M ,V ,PR ,VOL) 

REAL M,N,KI1 ,IT 

SF=4.0 

UT=850. 

PHIT=0.5 
VA=0 .6 
CL=2.0 
P 1=3.14159 
RH0M=15.12 
ALPHA=1 .2 
Al=l 2.0 
Xl= .05 
KI1=4. 93E-5 
THETA=0.8 
SIGMAY=12.54E6 

RT=SQRT(M*V/(PHIT*UT*PI* ( 1 . -V A**2) ) ) 
YB=SQRT(1.-PR**(-1./G)*(1.-VA**2) ) 

Cl= . 2* ( Al+1 . ) *1 . /V* ( 1 . -VA**2)*UT**2*RT**2*PHIT**2*X 1 /K 1 1 
C2=-l*THETA*RH0M*UT**2/2. 

C3=2.*SIGMAY/(SF*ALPHA) 

H=( VB-VA ) /20. 

SUM=0 .0 

X1=0.0 

Y1=VA 

DO 100 J=1 ,21 

N=C1*( (1.-Y1)/(1.+Y1) )*(PHIT**2 + .49*Y1**2)/(PHIT**2+.7225*Y1**2 
*) 

D=C3+C2*(1.-Y 1**2) 

DN=C1*(-1.*(PHIT**2+.49*Y1**2 ) / (PHIT**2+.7225*Y1**2 )*(1 .+Y1) - ( 1 
*.-Yl)*(l.+Yl)**(-2)*(PHIT**2+.49*Yl**2)/(PHIT**2+.7225*Y 1**2)+ 
*.98*Y1* (1 . -Yl) / (1 ,+Yl) / (PHIT**2 + .7225*Y1**2 )-l .445*Y 1* (1 . -Yl) / 
*(l.+Yl)*(PHIT**2+.49*Yl**2)*(PHIT**2+.7225*Yl**2)**(-2) ) 

DD=-2. *C2*Y1 

X2=.5*(1./(N/D)**.5)*( (DN*D-DD*N)/D**2) 

I F ( X 1 .EQ. 0.0) GO TO 50 
SUM=SUM+(Xl+X2)/2. 



6 9 



50 X1-X2 

Y1=Y1+H 
100 CONTINUE 
IT=-H*$UM 

V0L=2.*P I*RT**2*CL*IT 

RETURN 

END 

SUBROUTINE REGEN (T2,T3,T5,T6,P,D,MB,PR,DPPREG,DHDREG,V0L) 

REAL MB 

T=T5/2+T6/2 

P1=P*PR 

SIGMA=ASIGMA(T,P1) 

IF(T2 .EQ. T3) GO TO 99 
DTLM=T3-T2 

GAMMA-1. /(6803.*P1)**.41*(.0254*D)**1. 18/ (.56*DTLM)**1. 41 
DH1=CP2(T5)-CP2(T6) 

P 1=3.14159 

AREG-P 1/4 .*SIGMA*GAMMA*(MB*.454)*(DH1*2326. )**1 .41*(1 .+DHDREG) 
***2.41*( l.+PR** (+2)*DHDREG** (-3) )** .41*35 .3 
VOL=AREG*DPPREG** (- .41) 

99 RETURN 
END 

SUBROUTINE RECUP (T2H,T3H,T5H,T6H,PRH,DPPH,D,MA,MF,DPPREC,V0L) 

REAL MA,MF 

T=T5H/2.+T6H/2. 

PC=1./(1.-DPPH) 

P=14. 7*PRH 
SI GMA=ASIGMA(T,P) 

DTLM=T3H-T2H 

GAMMA- (6803. *P )**(-. 41) *( .0254*D)**1 . 18* ( .56*DTLM)**(-1. 41) 
DH1=CP2(T5H )-CP2(T6H) 

PI-3.14159 
DHDREC- .65 
CM-MA+MF 

AREC=P I/4.*SIGMA*GAMMA*(CM*.454)*(DH1*2326. )**1 .41* (1 ,+DHDREC) 
***2 .41* (l.+DHDREC**(-3)*PC**2)** .41*35 .3 
VOL=AREC*DPPREC**(- .41) 

RETURN 

END 

SUBROUTINE COOLER (T7,T1,P,D,NTP,MB,DHDCLR,DPPCLR,V0L) 

REAL NPP ,NTP,MB 

IF ( T7 .LT. Tl) GO TO 99 

TA-T7/2.+T1/2. 

SIGMA=ASIGMA(TA,P ) 

TC1=530. 

TC2-540. 

D TLM- ( Tl -T Cl -T7 +T C2 ) /L OG ( ( Tl -T Cl ) / ( T7 -T C2 ) ) 

GAMMA- (6803. *P)** (- .41)*( .0254*D)**1 ,18*( .56*DTLM)** (-1 .41 ) 
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DH1=CP2(T7)-CP2(T1 ) 

P 1=3. 14159 
TW=TCl/2.+TC2/2. 

RK=AK(TA,P)/l.E-4 
RMU=AMU(TA) /6.4E-4 
RM=(T7-T1)*ACP(TA)/(1 .*10. ) 

RPR=APR(TA,P)/6.4 
NPP= .80 

RRH0=ARH0(TA,P)/64.0 
RSV=ARH0(T7,P)/64.0 
ALPH A=RK*RMU** ( - . 8 ) *RM** . 8*RPR** . 4 
CB=RM**(-2.75)*RMU**(-.25)*RRH0*RSV/NTP/NPP 
ACLR=PI/4.*SIGMA*GAMMA*(MB* .454)*(DH1*2326. )**1 .41* (l.+DHDCLR) 
**(l.+ALPHA*DHDCLR)**1.41*(l.+CB*DHDCLR**(-3) )**.41*35.3 
V0L=ACLR*DP PCLR** (-.41) 

99 RETURN 
END 

SUBROUTINE COMBUSTOR (T3 ,T4,T4H ,T5H ,WP ,PR ,D ,MB ,PRH,MA,MF ,VOL) 
REAL MB,MA,MF,M 

COMMON /C0MB/DPPRG,DPPRC,DPPHE,DHDHR,DPHE2 

T=T3/2+T4/2 

P=WP*(PR-DPPRG) 

S I GM A=AS I GI'I A ( T , P ) 

DTLM=(T5H-T3-T4H+T4)/L0G( (T5H-T3) /(T4H-T4) ) 
GAMMA=(6803.*P)**(-.41)*( .0254*D)**1 .18*( .56*DTLM)** (-1 .41 ) 
DH1=CP2(T4)-CP2(T3) 

P 1=3 .14159 
TH=T4H/2+T5H/2 
TRH=(1. -DPPHE)*PRH 
PH=14.7* (PRH-DPPRC) 

RK=AK(T,P)/AK(TH,PH) 

RMU=AMU (T) /AMU (TH) 

M=MA+MF 

RM=MB/M 

RPR=APR(T,P)/APR(TH,PH) 

ALPHA=RK*RMU**(- .8) *RM**.8*RPR**.4 

AHTR=PI/4.*SIGMA*GAMMA*(MB*.454)*(DH1*2326. )**1 .41*( 1 .+DHDHR) * 
* (1 . +ALPHA*DHDHR)**1 .41*DHDHR**(-1.23)*35.3 

V0L=AHTR*DPHE2**(-.41) 

RETURN 

END 

SUBROUTINE WHBVOL (MB,MS,T6 ,T7 ,NWHB ,H2S ,H3S ,DPPWHB ,VOL) 

REAL MB-,MS,NWHB ,NTU 

T=T6/2+T7/2 

CP=ACP(T) 

CPS=2.0 

CR=MS*CPS/(MB*CP) 

I F ( CR .GT. 1.) GO TO 10 
CMIN=MS*CPS 
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GO TO 20 
10 CR=1./CR 
CMI N=CP*MB 
20 PS=4.0 

ER= ( NWHB*CR- 1 . ) / ( NWHB- 1 . ) 

NTU=(-PS)*L0G(1.+1./CR*L0G(1.-CR*(1.-ER**(1./PS))/(CR-ER**(1./ 

*PS)))) 

DH=H3S-H2S 

VOL= .00045*CMIN*NTU*MS*DH*DPPWHB**(- .41) 

RETURN 

END 

SUBROUTINE CONDVOL (MS, X, VOL) 

REAL MS 

V0L=MS*X*27.3 

RETURN 

END 
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SUBROUTINE STMTAB (NI.XI ,TI) 

DIMENSION SATL I Q ( 3 ) ,STEAM(3) , COEFFT (237 ) 

COMMON /STM/PRESS, TEMP ,QUAL, SATL IQ, STEAM 
DATA (COEFFT( I) ,1=1,30) / .22082221E+04, .65000000E+03 

.60464261E+01, .70000000E+01, .18552680E+00 
- .52279952E+00 



.69573391E+03 

- .19513682E+00 

- .50625849E+00 
.55592942E+03 
•26548914E-04 
.30706919E+00 
.45000000E+03 

DATA (COEFFT ( I ) 

. 32893178E- 1 4 
.81536841E-02 
.6427 2696E- 2 1 
.64251564E-10 
.33953897E+01 
.9275654 3E-17 
.39936598E-08 
DATA ( COEFFT ( I ) 
-.147 22587E- 2 1 
.27 45550 5E— 1 2 
.98512754E+00 
.73857 180E-22 
- .70679233E-14 
•30734203E-04 

- .33608032E-23 
DATA (COEFFT ( I ) 

. 1 5587770E-07 
.55000000E+03 
. 44979672E-1 2 
. 123996 14E-02 
. 1 294537 IE-21 
.1754800 IE-12 
.10150974E+01 
DATA ( COEFFT ( I) 
.6921785 5E — 1 6 
.59598410E-08 
.25000000E+03 
.49239634E-17 
.85546496E-06 
.19000000E+01 
•18796109E+03 
DATA (COEFFT ( I ) 
.11199490E+04 
.62379937E+05 
•35122466E+04 
.30000000E+01 
.52122736E+03 
.12340954E-04 
.18000000E+01 



.72595720E+00 
.20000000E+01 
.91349306E-04 
•36241343E+01 
.55000000E+03 
1= 31, 60) / 

-.117 34466E— 1 1 
•69519265E+01 
.12704366E-18 
.4479732 2E-07 
.43010153E+03 

- .95470656E-15 

- . 99567313E- 06 
1= 61, 90) / 

. 66553573E- 1 9 

- .18335302E-10 
.25000000E+03 
.3322550 6E-1 9 
.901 1420 IE- 1 1 
.21762691E-01 
•85819792E-20 

1= 91,120) / 

•81264440E-05 
•94743147E-18 
- . 51337380E-1 0 
.1279855 1E+01 
- . 6 5873139E- 1 9 
.23202424E-09 
.21848036E+03 
1=121,150) / 

- . 1 8345822E- 1 4 
.45575915E-06 
.50084539E-23 
.1625035 IE-13 
•14285709E-02 

- .1108507 7E+05 
.15699797E+03 

1=151,180) / 

.15000000E+01 
.29474498E+05 
-.12340833E+04 
•13000000E+04 
-.11 682947E- 01 
.13555136E-03 
. 2382209 1E+02 



.91 1489 58E-01 
.14585694E+02 

- .74133440E-06 
.1 5679922E-02 
.44960616E+02 

- .98736888E-22 
.30041 189E-1 7 
.1 1244476E-07 
.25000000E+03 

- .9697616 3E - 1 6 

- .28104756E-04 
.54927377E+03 

- .2900467 7E-1 2 
•59157992E-03 
.55000000E+03 
•47725605E-17 

- .31160333E-06 
.47000000E+03 
.3007557 5E-18 
•58149642E-09 
.25000000E+03 
.257 29365E- 1 8 
.60506309E-13 
. 1 7002639E-01 
•24971478E-16 
.49178986E-08 
.5492737 7E+0 3 

- . 1 3817 574E- 1 7 
•36145967E-06 
.45000000E+03 
.21134438E-20 

- .66417 458E- 1 2 
.12325099E-02 

-.15659070E-21 

- .74068320E-12 
.36753405E+00 

- .15994587E+05 

- .20213342E+03 

.15040455E+03 

- .14863585E+06 

- .18907586E+05 

- .68928827E+02 

.70848990E+01 

-.12001588E+00 

- .1 509931 9E-03 
.18776649E+03 



•61843682E-01 
.12404319E+03 
. 7332432 3E-06 
.25032530E-01 
.17926732E+03 
. 1 9301230E- 19/ 

- .421 13669E-1 5 

- .78419653E-05 

.20575548E-23 
•83370000E-13 
.1749054 1E-0 1 

- .18444779E-19 
.3535329 2E-1 0 
.78143274E+00/ 
.21848036E+03 

- .53349163E-14 

- .11956354E-03 
.55000000E+03 

- . 7 389865 IE- 1 5 
.867 12981E- 07 

- .78459129E-25 

- .3159 1904E- 1 5/ 
•13684599E-10 
.45000000E+03 

- .38302588E-14 
•48326895E-05 
•25000000E+03 
•62520106E-14 
•12506461E-03 
.55000000E+03/ 
.15327038E— 18 

. 14589563E- 10 
.74969797E+00 

- .460991 19E- 18 
•10479899E-08 
.17000000E+01 
.91407688 E+0 4 

-.48650176E+02/ 

- .19050547E+03 

-.13414819E+06 

- .5785887 6E+03 
.12036836E+04 
.16507788E+02 
.20411991E+01 
.30000000E+01 
.33080504E+03/ 
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DATA (COEFFT(I), 1=181, 210) /- .32141408E+01, .42449101E+03 

*, .16400835E+04, .18746657E+02, .63978878E+02 , .20719841E+04 
*, .13000000E+04, .18000000E+01, -.64369800E-05, .96418803E-02 
*, .38879248E+01 , .11034060E-05, .47225001E-03 .94113600E+01 
*, .69416973E-05,-.60892120E-02, .10951922E+01, .13000000E+04 
*, .18000000E+01 .33101208E-03 , .21661916E+01 , .53500291E+03 
.77236618E-03, .36753650E+00,-.25450121E+03, .25249707E-02 
.16427106E+01 , .61477766E+03, 0.0, 0.0 / 

DATA (COEFFT(I), 1=211, 237) / 0.0, 0.0, 0.0, 0.0, 0.0 

*, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, .9000000QE+03 
*, .14000000E+01,-.15478259E-04, .19351556E-01, .24698909E+01 
*, .23329644E-04, .1149379 IE-02 , - .13139208E+02 , .58860219E-05 
*,-.12592899E-01, .55101160E+01/ 

N=NI 

X=XI 

Y=Y I 

I=N-3 

LINK=0 

PRESS=X 

DELTA1=0 .0 

ERR0R1=0.0 

ERR0R2=0.0 

IF (N-6) 1000,1000,1030 
1000 IF (I) 1020,1020,1010 
1010 L INK=1 
N=I 

1020 IF (X-COEFFT(N) ) 1040,1040,1620 
1030 N=I 

1040 GO TO (1130, 1060, 1190, 1230, 1080, 1090, 1180, 1330), N 
1050 X=Z 
N=2 

1060 TEMP=X 

GO TO 1190 
1070 PRESS=EXP(Z) 

GO TO 1150 
1080 11=2 
JJ=3 

TOLER= .005 
GO TO 1100 
1090 11=3 
JJ=2 

TOLER= .000005 
1100 H=Y 
1110 LINK=2 
1120 N=1 

1130 X=ALOG(PRESS) 

PLOG=X 
GO TO 1190 
1140 TEMP=Z 
1150 X=TEMP 
J=1 
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N=4 

GO TO 1190 
1160 SATLIQ(J)=Z 
J=J+1 
N=N+1 

IF (J-3) 1190,1190,1170 
1170 IF (LINK) 1640,1640,1300 
1180 H=X 
S=Y 
X=S 

1190 I=23*N-19 
K=1 

IF (X-COEFFT(I)) 1200,1200,1210 
1200 K=12 
1210 I=I+K 

ARG=X-COEFFT( I ) 

Z=0.0 

DO 1220 K=1 ,10 
1 = 1+1 

1220 Z=Z*ARG+COEFFT ( I ) 

GO TO (1140,1070,1050,1160,1160,1160,1410) ,N 
1230 TEMP=Y 

GO TO 1300 

1240 ERROR=H-STEAM( 1 1 ) 

IF (ABS( ERROR) -TOLER) 1640,1250,1250 
1250 IF (ERR0R1) 1260,1280,1260 
1260 ERR0R1=ERR0R1 -ERROR 

IF (ABS(ERRORl)-TOLER) 1290,1290,1270 
1270 SL0PE=DELTA1/ERR0R1 
1280 DELTA1=SL0PE*ERR0R 
ERR0R1=ERR0R 
1290 TEMP=TEMP+DELT A1 
1300 IF (PRESS-COEFFT(l) ) 1320,1320,1310 

1310 IF ( ( .00001178*PRESS+.09411)*PRESS+382. 1-TEMP) 1320,1320,1620 
1320 T=.55555556*TEMP+255. 38223 
TLOG=ALOG(T) 

TAU=1.0/T 

COEFFT (222)=TAU*TAU*186210 .06 
C0EFFT(217)=EXP(C0EFFT(222)+7.8791476-TL0G) 

COEFFT (221X62460. *TAU 
C0EFFT(220)=126970.*TAU 
C0EFFT(219)=EXP (149. 27765-23. *TLOG) 

C0EFFT(212) =1. 89-C0EFFT(217) 

COEFFT (21 1 )=82 .546-COEFFT (22 1 ) 
C0EFFT(210)=.21828*T-C0EFFT(220) 

COEFFT (209 XOEFFT (219)- .000363 5*T 
C0EFFT(222)=(C0EFFT(222)+C0EFFT(222)+1.0) *C0EFFT(2 17) 

COEFFT (217 )=COEFFT (212 )-COEFFT (222) 

C0EFFT(2 14) =C0EFFT(217) /C0EFFT(2 12) 

COEFFT (216)=C0EFFT(211)*C0EFFT (214)+41 .273-COEFFT (221) 
C0EFFT(215)=C0EFFT(210)*C0EFFT(214)-.5*C0EFFT(220) 
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COEFFT (214) =C0EFFT(209)*C0EFFT(2 14)+1 .8461538*C0EFFT(2 19) 

COEFFT (221)= .50*C0EFFT (211 )-COEFFT (216) 
C0EFFT(220)=.25*C0EFFT(210)-C0EFFT(2 15) 

COEFFT (219)= .O76923O77*C0EFFT (209 )-COEFFT (214) 

C0EFFT(218)=( ,0037299965*T+14. 531813) *T+1 7806. 5 13+472. 2493*TL0G 
C0EFFT(223)=-.0074599931*T-14.531813*TL0G+2. 4524439+472. 24937 
**TAU 

1330 P=.068046190*PRESS 
PLOG=ALOG(PRESS) 

COEFFT (213)=C0EFFT(212) *TAU 

C0EFFT(226)=C0EFFT(2 13) *P 

COEFFT (225) =COEFFT (226) *COEFFT (226) 

C0EFFT(224) =C0EFFT(225) *C0EFFT(225) 

COEFFT (224)=C0EFFT (224)*C0EFFT (224) *COEFFT (226 ) 

C0EFFT(226) =C0EFFT(226)*C0EFFT(213) 

COEFFT (213) =P*TAU 

STEAM(1 )=( ( (C0EFFT(209)*C0EFFT(224)+C0EFFT(210) )*C0EFFT(225)+ 
*COEFFT (211 ) )*COEFFT (226 ) +COEFFT (212 )+4 .55504/C0EFFT (213) ) 
**.0160185 

STEAM(2 )=( ( ( (COEFFT(214)*COEFFT(224)+COEFFT(215) )*C0EFFT(225)+ 
*C0EFFT(216) )*COEFFT (226)+C0EFFT (217) )*P+COEFFT (218) )*. 04355685 
STEAM(3 )=( ( ( (C0EFFT(219)*C0EFFT(224)+C0EFFT(220) )*C0EFFT(225)+ 
*COEFFT (221) )*COEFFT (226)+C0EFFT (222) )*COEFFT (213) +C0EFFT (223) 
*+4.55504*PL0G)*(-. 02419825) 

QUAL=1 .0 
I=L INK+1 

GO TO (1640,1640,1340,1240,1340,1450) , I 
1340 X=STEAM( 1 1 ) 

IF (H-X) 1350,1640,1360 
1350 Y=SATLIQ( 1 1 ) 

QUAL=(H-Y ) / (X -Y ) 

X=1 .O-QUAL 

STEAM(1 )=QUAL*STEAM(1 )+X*SATLIQ(l ) 

STEAM (J J ) =Q UAL*STE AM ( J J ) +X*SATL I Q ( J J ) 

STEAM( II) =H 

IF (LINK-2) 1640,1640,1460 
1360 1= ( 1 1-1 ) *1 1+154 
J=1 

X=PLOG 

1370 DO 1380 K=209 ,219 

COEFFT(K)=COEFFT( I ) 

1380 1=1+1 

X=X-C0EFFT(209) 

Y=Y-COEFFT (210) 

C0EFFT(209)=(C0EFFT(211)*X+C0EFFT(212) )*X+C0EFFT(213) 

COEFFT (210 )=( COEFFT (21 4) *X+COEFFT (215)) *X+COEFFT (216) 

COEFFT ( 2 1 1 ) = ( (C0EFFT(217) *X+C0EFFT(2 18) )*X+C0EFFT(219) )*Y 
Z=COEFFT (211 )+COEFFT (210) 

SLOPE=COEFFT(211)+Z 
Z=Z*Y+COEFFT (209) 

GO TO (1390,1430,1530,1540) ,J 
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1390 LINK=3 

IF (Z-TEMP) 1300,1300,1400 
1400 TEMP=Z 

GO TO 1300 
1410 J=2 
1=227 
X=h 

IF (Z-H) 1420,1370,1370 
1420 1=187 

GO TO 1370 
1430 DPDS=SLOPE 
PRESS=EXP( Z) 

PLOG=Z 

IF (1-227) 1520,1520,1440 
1440 LIMK=4 
11=2 
JJ=3 

GO TO 1120 
1450 DH=H-STEAM(2) 

STEAM(3)=STEAM(3)+DH/ (TEMP+459.69) 

1460 DS=S-STEAM(3) 

DELTA1=STEAM(3)-ERR0R1 
IF (ERR0R1) 1470,1490,1470 
1470 IF (ABS(DELTAl)-. 000005) 1500,1500,1480 
1480 DPDS=(PL0G-PL0G1) /DELTA1 
1490 PL0G1=PL0G 

ERR0R1=STEAM(3) 

1500 PRESS=PRESS*(1.0+DS*DPDS) 

IF (LINK-4) 1510,1510,1550 
1510 IF (ABS(DS)-. 000005) 1640,1120,1120 
1520 LINK=5 
X=H 
Y=S 
1=198 
J=3 

GO TO 1370 
1530 TEMP=Z 

DTDS=SLOPE 

X=PLOG 

Y=H 

1=165 

J=4 

GO TO 1370 
1540 DTDH=SLOPE 
GO TO 1300 

1550 TEMP=TEMP+DS*DTD$ 

IF (ABS(DH )- .005) 1610,1560,1560 
1560 DELTA2=DH-ERR0R2 

IF (ERR0R2) 1570,1590,1570 
1570 IF ( ABS(DELTA2 )- .005) 1600,1600,1580 
1580 DTDH=(TEMP1-TEMP)/DELTA2 
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1590 TEMP1=TEMP 
ERR0R2=DH 

1600 TEMP=TEMP+DH*DTD H 
GO TO 1320 

1610 IF (ABS(DS)-. 000005) 1640,1600,1600 
1620 WRITE (6,1630) NI,XI,YI 
1630 FORMAT ( 13 ,2E14 .8 ,13H OUT OF RANGE) 
STOP 

1640 RETURN 
END 
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APPENDIX B 
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BRAYTON CYCLE: PRESSURE RATIO =2. 00 

WORKING PRESSURE = 1 35 . 0 
TOTAL PRESSURE DR0P=.10 
MASS FLOW (Lbm/s) =444 . 37 



STATE 

1 

2 

3 

4 

5 

6 

PP 

7 

COMPONENT 

REGEN 

TURBINE 

COOLER 

COMPRESSOR 

HEATER 

BOILER 



TEMP (deg R) 

540.0 
675.7 

1668.8 

2060.0 

1804.2 

811.2 

784.0 
778.3 

DPP DHD 
. 078 . 744 

.007 .260 

.0005 

.014 



SUMMARY: THERMAL EFF=0.370 

HEATER EFF=0 .933 (PHI=0.90) 

SFC ( Lbm/HP-hr ) =0 . 572 

VOL TOTAL (-CU ft)=5l6.5 

POWER (HP) =25000.0 

WORK FRAC BY STM CYCLE=0.045 



VOL( cu ft) 

15.5 EFF=0 . 88 
10.0 EFF(P)=0.91 
139.8 

29.3 EFF( P) =0 . 88 
EFF=0 . 94 



HEATER CYCLE: PRESSURE RATIO =4.79 

TOTAL PRESSURE DR0P=.43 
AIR FLOW ( Lbm/s ) =47 . 84 
FUEL FLOW( Lbm/s )= 3*97 



STATE 


TEMP(deg R) 


lh 


535.0 


2h 


835.2 


3h 


1727.4 


4h 


5054.5 


5h 


1768.8 


6h 


945.0 


7h 


760.0 


COMPONENT 


DPP 


DHD 


RECUP 

TURBINE 


.321 


.650 


HEATER 

COMPRESSOR 


.012 


.310 



VOL( cu ft) 

2.6 EFF=0.96 

1.1 
6 7.4 

16.0 EFF( P) =0 . 85 



RANKINE CYCLE: SATURATION PRESSURE= 56.7 

STEAM FLOW ( Lbm/s )= 2.96 

FLOW RAT 10 =0.006 7 



STATE 


TEMP( deg R ) 


PRESSURE ( 


Is 


540.0 


0.5 


2s 


555-0 


56.7 


SF 


749.0 


56.7 


3s 


776.2 


56.7 


4s 


540.0 


0.5 


COMPONENT 


VOL( cu ft) 




BOILER 


162.3 




TURBINE 


4.3 EFF( IS ) =0 . 88 


CONDENSER 


68.3 





QUAL=0 . 844 
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BRAYTON CYCLE! PRESSURE RATIO = 3-00 
WORKING PRESSURE =135.0 
TOTAL PRESSURE DR0P=.10 
MASS FLOW (Lbm/s) =287 .63 



STATE 


TEMP( deg R) 


SUMMARY: THERMAL EFF=0.432 


1 


540.0 


HEATER EFF=0.933 (PHI=0.90) 


2 


769.8 


SFC(Lbm/HP-hr) =0.490 


3 


1538.8 


VOL TOTAL( cu ft)= 318.1 


4 


2060.0 


POWER( HP) =25000.0 


5 


1643.6 


WORK FRAC BY STM CYCLE=0.044 


6 


874.7 




PP 


841.2 




7 


831 .8 




COMPONENT 


DPP DHD 


VOL( cu ft) 


REGEN 


.050 .837 


11.5 EFF= 0 . 88 


TURBINE 




5.3 EFF(P) =0 .91 


COOLER 


. 024 . 352 


88.9 


COMPRESSOR 




20.9 EFF( P) =0 . 88 


HEATER 


.0004 




BOILER 


.026 


EFF=0 .94 



HEATER CYCLE: PRESSURE RATIO = 4.79 
TOTAL PRESSURE DR0P=.43 
AIR FLOW ( Lbm/s ) = 40.99 
FUEL FLOW ( Lbm/s ) = 3.4o 



STATE 


TEMP( deg R) 




lh 


535.0 




2 h 


835.2 




3h 


1586.5 




4h 


4932.9 




5h 


1638.8 




6 h 


945.0 




7h 


760.0 




COMPONENT 


DPP DHD 


VOL( cu ft) 


RECUP 


.310 .650 


2.2 EFF=0 .93 


TURBINE 




0.9 EFF( P ) =0 


HEATER 


.022 .358 


28.6 


COMPRESSOR 


12.7 EFF( P) =0 


RANKINE CYCLE: SATURATION 


PRESSURE= 128.0 


STEAM FLOW(Lbm/s) = 2.48 
FLOW RATI 0 = 0.0086 




STATE 


TEMP(deg R) 


PRESSURE (psi) 


Is 


540.0 


0.5 


2 s 


555.0 


128.0 


SF 


806.2 


128.0 


3s 


839.7 


128.0 


4s 


540.0 


0.5 QUAL= 


COMPONENT 


VOL(cu ft) 




BOILER 


88.3 




TURBINE 


3.4 EFF( IS ) =0 . 88 


CONDENSER 


55.5 
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BRAYTON CYCLE: PRESSURE RATIO = 4.00 

WORKING PRESSURE = 135-0 
TOTAL PRESSURE DR0P=.10 
MASS FLOW (Lbm/s) =2 39 . 91 



STATE 


TEMP( deg R) 


SUMMARY : THERMAL EFF=0.447 


1 


540.0 


HEATER EFF=0 .933 (PHI 


2 


843.9 


SFC(Lbm/HP-hr)=0.473 


3 


1453.8 


VOL TOTAL( cu ft)= 352.: 


4 


2060.0 


PO T WER( HP) =25000.0 


5 

6 

PP 

7 


1536.9 

927 .O 

876.8 

860.6 


WORK FRAC BY STM CYCLE 


COMPONENT 


DPP DHD 


VOL(cu ft) 


REGEN 

TURBINE 


.021 .930 


13.2 EFF=0 . 88 
3.6 EFF( P ) =0 . 91 


COOLER 

COMPRESSOR 

HEATER 


.032 .444 

.0003 


89.8 

18.4 EFF(P ) =0 . 88 


BOILER 


.046 


EFF=0 .94 


HEATER CYCLE: PRESSURE RATIO = 4.79 

TOTAL PRESSURE DR0P=.43 
AIR FLOW ( Lbm/ s ) = 39 • 59 
FUEL FLOW ( Lbm/s ) = 3.29 

STATE TEMP (deg R) 

lh 535-0 

2h 835.2 

3h 1494.5 

4h 4855-0 

5h 1553-8 

6 h 945 . 0 

7 h 760.0 


COMPONENT 


DPP DHD 


VOL(cu ft) 


RECUP 

TURBINE 


.299 .650 


2.1 EFF=0 . 92 
0.8 EFF(P)=0.89 


HEATER 

COMPRESSOR 


.033 .406 


17.7 

12.0 EFF( P) =0.85 


RANKINE CYCLE: SATURATION 

STEAM FLOW ( Lbm/s )= 3 . 19 


PRESSURE 1 199-9 


FLOW RAT 10=0. 


0133 




STATE 


TEMP(deg R) 


PRESSURE (psi) 


Is 


540.0 


0.5 


2s 


555.0 


199.9 


SF 


841.8 


199.9 


3s 


892.0 


199.9 


4s 

COMPONENT 

BOILER 


540.0 

VOL( cu ft) 
119.0 


0.5 QUAL=0 . 81 0 


TURBINE 

CONDENSER 


5.1 EFF( IS ) =0 . 88 
70.5 
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BRAYTON CYCLE: PRESSURE RAT 10= 5*00 

WORKING PRES SURE =135.0 
TOTAL PRESSURE DR0P=.10 
MASS FLOW (Lbm/s)=212.34 



STATE 


TEMP (deg R) 


SUMMARY: THERMAL EFF=0.460 


1 


540.0 




HEATER EFF=0 . 933 (PHI=0.90) 


2 


905.8 




SFC(Lbm/HP-hr)= 0.460 


3 


1391 .8 




VOL TOTAL( cu ft)=596.1 


4 


2060.0 




POWER (HP) =25000 . 0 


5 


1458.1 




WORK FRAC BY STM CYCLE= 0.111 


6 


972.1 






PP 


876.8 






7 


846.9 






COMPONENT 


DPP DHD 


VOL( cu ft) 


REGEN 


.010 .950 


15.7 


EFF=0 . 88 


TURBINE 




2.7 


EFF( P) =0 .91 


COOLER 


.022 .536 


107.8 




COMPRESSOR 




16.7 


EFF( P) =0 . 88 


HEATER 


.0003 






BOILER 


.068 




EFF=0 . 94 



HEATER CYCLE: PRESSURE RATIO= 4.79 

TOTAL PRESSURE DR0?=.43 
AIR FLOW ( Lbm/ s ) = 38.49 
FUEL FLOW ( Lbm/ s ) = 3*19 



STATE 


TEMP(deg R) 




lh 


535.0 




2h 


835.2 




3h 


1427.4 




4h 


4798.9 




5h 


1491.8 




6h 


945.0 




7h 


760.0 




COMPONENT 


DPP DHD 


VOL( cu ft) 


RECUP 


.289 .650 


2.1 EFF=0 .90 


TURBINE 




0.8 EFF( P) = 0 . 89 


HEATER 


.044 .454 


12.4 


COMPRESSOR 




11.5 EFF(P)=0 .85 


RANKINE CYCLE: SATURATION 

STEAM FLOW(Lbm/s) = 5-21 


PRESSURE = 200.0 


FLOW RATI0=0 


.0245 




STATE 


TEMP (deg R) 


PRESSURE(psi) 


Is 


540.0 


0.5 


2s 


555.0 


200.0 


SF 


841.8 


200.0 


3s 


937.1 


200.0 


4s 


540.0 


0.5 QUAL=0 . 826 


COMPONENT 


VOL( cu ft) 




BOILER 


298.1 




TURBINE 


10.9 EFF( IS ) =0. 88 


CONDENSER 


117.4 
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BRAYTON CYCLE: PRESSURE RATIO = 6.00 

WORKING PRESSURE =135.0 



TOTAL PRESSURE DR0P=.10 
MASS FLOW(Lbm/s) =197 .00 




STATE 


TEMP( deg R) 


SUMMARY: THERMAL EFF= 0.464 


1 


540.0 


HEATER EFF=0 .933 (PHI 


2 


959.4 


SFC(Lbm/HP-hr) =0.456 


3 


1343.8 


VOL TOTAL( cu ft) =879-0 


4 


2060.0 


POWER ( HP) =25000.0 


5 


1396.2 


WORK FRAC BY STM CYCLE 


6 


1011.8 




PP 


876.8 




7 


8350 




COMPONENT 


DPP DHD 


VOL( cu ft) 


REGEN 


.010 .950 


14.4 EFF=0 .88 


TURBINE 




2.2 EFF ( P ) = 0 . 9 1 


COOLER 


.016 .628 


129.7 


COMPRESSOR 




15.8 EFF( P) =0 . 88 


HEATER 


.0003 




BOILER 


.074 


EFF=0 .94 


HEATER CYCLE! 


: PRESSURE RATIO= 4.79 


TOTAL PRESSURE DR0P=.43 
AIR FLOW ( Lbm/ s ) = 38 . 1 7 




FUEL FLOW ( Lbm/s ) - 3 . 17 




STATE 


TEMP( deg R) 




Ih 


535.0 




2h 


835.2 




3h 


1375.4 




4h 


4755.8 




5h 


1443.8 




6h 


945.0 




7h 


760.0 




COMPONENT 


DPP DHD 


VOL( cu ft) 


RECUP 


.278 .650 


2.1 EFF= 0 . 89 


TURBINE 




0.8 EFF ( P ) = 0 . 89 


HEATER 


.055 .502 


9.5 


COMPRESSOR 




11.4 EFF( P) =0 . 85 


RANKINE CYCLE : SATURATION 

STEAM FLOW (Lbm/ s ) = 6.70 


PRESSURE= 200.0 


FLOW RATIO =0. 


.0340 




STATE 


TEMP( deg R) 


PRESSURE(psi) 


Is 


540.0 


0.5 


2s 


555.0 


200.0 


SF 


841.8 


200.0 


3s 


976.8 


200.0 


4s 


540.0 


0.5 QUAL=0 .838 


COMPONENT 


VOL(cu ft) 




BOILER 


523.6 




TURBINE 


16.2 EFF( IS ) =0 . 88 


CONDENSER 


153.3 
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BRAYTON CYCLES PRESSURE RATI0=7.00 
WORKING PRESSURES 35.0 
TOTAL PRESSURE DR0P=.10 
MASS Flow(Lbm/s)=187.48 



STATE 


TEMP (deg R) 


SUMMARY; THERMAL EFF=0.463 


1 


540.0 


HEATER EFF=0 .933 (PHI=0.90) 


2 


1007.0 


SFC (Lbm/HP=hr ) =0 .457 


3 


1304.9 


VOL TOTAL( cu ft) =1195. 7 


4 


2060.0 


POWER (HP) =25000.0 


5 


1345.5 


WORK FRAC BY STM CYCLE =0.1 89 


6 


1047.6 




PP 


876.8 




7 


825.2 


VOL( cu ft) 


COMPONENT 


DPP DHD 


REGEN 


.010 .950 


13.6 EFF=0 . 88 


TURBINE 




1.9 EFF( P) =0 .91 


COOLER 


.012 .720 


155.9 


COMPRESSOR 




15.3 EFF( P) =0 . 88 


HEATER 


.0004 




BOILER 


.078 


EFF=0 .94 


HEATER CYCLE; 


PRESSURE RATIO = 4.79 


TOTAL PRESSURE DR0P=.43 




AIR FLOW(Lbm/s) = 38.22 
FUEL FLOW(Lbm/s)= 3-17 




STATE 


TEMP( deg R) 




lh 


535.0 




2h 


835.2 




3h 


13330 




4h 


4721.2 




5h 


1404.9 




6h 


945.0 




7h 


760.0 




COMPONENT 


DPP DHD 


VOL( cu ft) 


RECUP 


.267 .650 


2.1 EFF=0 . 87 


TURBINE 




0.8 EFF ( P) =0 . 89 


HEATER 


.066 .550 


7.7 


COMPRESSOR 




11.4 EFF( P) =0 . 85 


RANKINE CYCLE: SATURATION 


PRESS URE= 200.0 


STEAM FLOW ( Lbm/ s ) = 7. 92 
FLOW RATI0=0.0422 




STATE 


TEMP( deg R) 


PRESSURE (psi) 


Is 


540.0 


0.5 


2s 


555.0 


200.0 


SF 


841.8 


200.0 


3s 


1012.6 


200.0 


4s 


540.0 


0.5 QUAL=0 . 849 


COMPONENT 


VOL( cu ft) 




BOILER 


782.4 





TURBINE 21.0 EFF(IS)=0.88 

CONDENSER 18 3. 6 
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BRAYTON CYCLE 


I: PRESSURE RAT 10= 8.00 


WORKING PRESSURE = 135.0 




TOTAL PRESSURE DROP= .10 




MASS FL0W(Lbm/s) = 181.22 




STATE 


TEMP (deg R) 


SUMMARY: THERMAL EFF=0.460 


1 


540.0 


HEATER EFF=0.933 (PHI 


2 


1049.8 


SFC(Lbm/HP-hr) =0.459 


3 


1272.5 


VOL TOTAL( cu ft) =1550 


4 


2060.0 


POWER (HP) =25 000.0 


5 


1302.9 


WORK FRAC BY STM CYCL! 


6 


1080.2 




PP 


876.8 




7 


816.3 




COMPONENT 


DPP DHD 


VOL( cu ft) 


REGEN 


.010 .950 


13.1 EFF=0 . 88 


TURBINE 




1.6 EFF( P) =0 .91 


COOLER 


.010 .812 


181 .1 


COMPRESSOR 




15.1 EFF( P) =0 . 88 


HEATER 


.0004 




BOILER 


.080 


EFF=0 .94 


HEATER CYCLE: 


PRESSURE RAT 10= 4.79 


TOTAL PRESSURE DR0P=.43 




AIR FLOW (Lbm/s) =38.45 




FUEL FLOW ( Lbm/s ) = 3 . 19 




STATE 


TEMP (deg R) 




lh 


535.0 




2h 


835.2 




3h 


1298.2 




4h 


4692.5 




5h 


1372.5 




6h 


945.0 




7h 


760.0 




COMPONENT 


DPP DHD 


VOL ( c u ft) 


RE CUP 


.256 .650 


2.1 EFF=0 . 86 


TURBINE 




0.8 EFF( P) = 0 . 89 


HEATER 


.076 .598 


6.5 


COMPRESSOR 


11.5 EFF(P)=0 .85 


RANKINE CYCLE 


: SATURATION 


PRESSURE=200 . 0 


STEAM FLOW ( Lbm/s )= 8.97 




FLOW RATIO =0. 


0495 




STATE 


TEMP( deg R) 


PRESSURE (psi) 


Is 


540.0 


0.5 


2s 


555-0 


200.0 


SF 


841.8 


200.0 


3s 


1045.2 


200.0 


4s 


540.0 


0.5 QUAL=0 . 859 


COMPONENT 


VOL( cu ft) 




BOILER 


IO 83 .I 




TURBINE 


25.5 EFF( IS ) =0. 88 


CONDENSER 


210.2 
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BRAYTON CYCLE: 


PRESSURE RATIO= 9.00 


WORKING PRESSURE =135.0 




TOTAL PRESSURE 


DROP= . 1 0 




MASS FLOW (Lbm/s)=l 76 . 95 




STATE TEMP( deg R) 


SUMMARY: THERMAL EFF=0.456 


1 


540.0 


HEATER EFF=0.933 


2 


1089.0 


SFC(Lbm/HP-hr)=0.‘ 


3 


1245.0 


VOL TOTAL( cu ft)= 


4 


2060.0 


POWER ( HP) =25000.0 


5 


1266.2 


WORK FRAC BY STM 1 


6 


1110.2 




PP 


876.8 




7 


808.2 




COMPONENT 


DPP DHD 


VOL( cu ft) 


REGEN 


.010 .950 


12.7 EFF= 0 . 88 


TURBINE 




1.4 EFF ( P ) = 0 . 9 1 


COOLER 


.010 .904 


197.5 


COMPRESSOR 




15.0 EFF( P) = 0 . 88 


HEATER 


.0004 




BOILER 


.080 


EFF= 0 .94 


HEATER CYCLE: 


PRESSURE RATIO= 4.79 


TOTAL PRESSURE 


DROP = .43 




AIR FLOW(Lbm/s) 


1 = 38 . 8 O 




FUEL FLOW ( Lbm/ s ) = 3.22 




STATE TEMP (deg R) 




lh 


535.0 




2h 


835.2 




3h 


1268.3 




4h 


4668.2 




5b 


1345.0 




6h 


945.0 




7h 


760.0 




COMPONENT 


DPP DHD 


VOL( cu ft) 


RECUP 


.245 .650 


2.2 EFF=0 . 85 


TURBINE 




0.8 EFF( P) =0 . 89 


HEATER 


.087 .646 


5.7 


COMPRESSOR 




11.7 EFF(P)=0.85 


RANKINE CYCLE: 


SATURATION 


PRESSURE 3 200.00 


STEAM FLOW ( Lbm/s ) = 9-90 




FLOW RATIO= 0.0560 




STATE TEMP( deg R) 


PRESSURE (psi) 


Is 


540.0 


0.5 


2s 


555.0 


200.0 


SF 


841.8 


200.0 


3s 


1075.2 


200.0 


4s 


540.0 


0.5 QUAL=0.867 


COMPONENT 


VOL( cu ft) 




BOILER 


1439 .0 




TURBINE 


29.8 EFF(IS) =0.88 


CONDENSER 


234.5 
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BRAYTON CYCLE: PRESSURE RATI0=10,00 

WORKING PRESSURES 35.0 
TOTAL PRESSURE DR0P=.10 
MASS FLOW (Lbm/s)=174.00 



STATE 

1 

2 

i 

6 

PP 

7 



TEMP(deg R) 
54o.o 

1125.1 

1221.1 

2060.0 

1234.2 

1138.1 
876.7 
800.9 



SUMMARY: THERMAL EFF=0.451 

HEATER EFF=0 .933 (PHI=0.90) 

SFC ( Lbm/HP-hr ) =0 . 469 

VOL TOTAL ( cu ft) =2403.0 

POWER (HP) =25000.0 

WORK FRAC BY STM CYCLE=0.296 



COMPONENT 


DPP 


DHD 


VOL(cu ft) 


REGEN 


.010 


.950 


12.4 


EFF=0 . 88 


TURBINE 






1.3 


EFF( P) =0 .91 


COOLER 


.008 


.996 


232.8 




COMPRESSOR 






14.9 


EFF( P) =0 . 88 


HEATER 


.0004 








BOILER 


.082 






EFF=0 .94 



HEATER CYCLE: PRESSURE RATIO= 4.79 

TOTAL PRESSURE DR0P=.43 
AIR FLOW ( Lbm/ s ) - 39-22 
FUEL FLOW(Lbm/s) = 3.25 



STATE 


TEMP (deg R) 






lh 


535-0 






2h 


835.2 






3h 


1242.5 






4h 


4647.3 






5h 


1321 .1 






6h 


945.0 






7h 


760.0 






COMPONENT 


DPP DHD 


VOL( cu : 


ft) 


RECUP 


.235 .650 


2.2 


EFF=0 . 84 


TURBINE 




0.8 


EFF( P) =0 . 89 


HEATER 


.098 .694 


5-0 




COMPRESSOR 




11 .9 


EFF(P) =0 .85 



RANKINE CYCLE: SATURATION PRESSURE= 199.9 

STEAM FLOW ( Lbm/ s ) = 10 . 76 
FLOW RATI0=0 . 06 l 8 



STATE 


TEMP( deg R) 


Is 


540.0 


2 s 


555-0 


SF 


841.7 


3s 


II 03 .I 


4s 


540.0 


COMPONENT 


VOL( cu ft 


BOILER 


1830.6 


TURBINE 


34.0 e: 


CONDENSER 


257.1 



PRESSURE (psi) 

0.5 

199.9 

199.9 

199.9 

0.5 QUAL=0.875 



IS ) =0 . 88 



88 
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